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PREFACE 

From time to time scientific and engineering personnel find themselves in 
situations where there is a need to know how to "program" their problems for the 
modern electronic digital computer. 

To learn machine language or assembly language programming requires too much 
time; there are too many petty rules, and numerous exceptions are involved. Unless 
the scientist can devote full time to the effort, he looks down the “long road ahead® 
and is discouraged from even trying to start. 

This short course is an attempt to give to the scientist and engineer, sufficient 
skill, in a minimum of time, to enable him to efficiently program his own problems. 
To accomplish this objective, the course contains the following characteristics: 

~ « The FORTRAN language is used exclusively! 
- Non essential concepts have been eliminated! 


- “Learning Machine" techniques are used--in that continued lesson successes 
depend upon learning previous concepts! 


. Every effort (through “humor", etc) is made to keep you awake! 


GETTING STARTED? 

eo 0 6 © © © 0 © co ce co oe 18 SO easy? 
SETTLE DOWN? 

ee © © © eo © oe oe o in anything? 
SCRATCH PAPER? 


° oe 0 ee wo oo © oe YOS, have some at hand! 


FORTRAN MANUALS? 


© 0 0 eo 0 © © © ow oe o NOt needed, but available! 


IS THAT ALL? 


oe ee ee eo 8 wo oo Yed, go to.Card 1, 


“Her name is Constance", said the man, “my ever faithful and unchanging wife". 
In a similar vein, the mathematician likes to call unchanging magnitudes, CONSTANTS. 
How many constants would you say there are among the following? 


a) 1062 

b) 3/4 

c) 1062. 

d) .52E2 ("E2" indicates 2 is the exponent of 10) 


If your answer is 4, go to 3 
If your answer is less than 4, go to 2 


You probably thought the fourth example was not a CONSTANT because of the 
letter "E”, This is a good guess, but it is too cautious! The "E* is just a symbol 
to indicate "exponent of 10". Some constants are more easily expressed this way. 
For example, the constant. 


62,000, 000,000,000,000,000,000,000,000,000 = .62E32 
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Now that your concept of CONSTANTS is 0. K., do you see a basic format 
difference in the constants, 1062 and 1062.7? Yes, there is a decimal point expressed 
in one. Fortran makes use of two types of constants, FIXED POINT CONSTANTS and 
FLOATING POINT CONSTANTS. 


Floating point format makes it possible to do arithmetic (add, subtract, etc.) 
and have the computer automatically position the decimal point, thus making it 
unnecessary for you to continually scale and rescale as the problem progresses. The 
compiler must know what type of CONSTANT you are thinking about. Therefore. 

IF THE DECIMAL POINT IS EXPRESSED, YOU ARE USING FLOATING POINT CONSTANTS. 


IF THERE IS NO DECIMAL POINT EXPRESSED, YOU ARE USING FIXED POINT CONSTANTS, 


Go To & 


A plus (+) or minus (-) sign preceding a constant indicates a positive or 
negative constant. A plus (+) need not be written with positive constants, but the 
minus (-) must be written with negative constants. In summary then, if you use. 


-6231E-7 
etc. 
you are thinking 


about FLOATING 
POINT CONSTANTS! 


+2 
61235 
etc. 
you are thinking 


about FIXED 
POINT CONSTANTS?, 
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35 
23 
35 


you are thinking 
about "LIVE” 
CONSTANCES ! 


In addition to constants, you will make use of variables. VARIABLES are 
quantities that may assume a succession of changing magnitudes. Whereas the 
representation of each constant is the constant itself, each variable mst be 
represented by a combination of letters am digits. 


FIXED POINT VARIABLES only represent integers or whole numbers--expressed without a 
decimal point 


FLOATING POINT VARIABLES represent fractional, mixed,or whole numbers with a decimal 
point expressed. 
The following fqur rules apply for variables: 


(1) The first character of the variable name mst be alphabetic. 


(2) A first character of I, J, K, L, M, or N indicates a FIXED POINT VARIABLE; any 
other first letter indicates a FLOATING POINT VARIABLE, 


(3) Each variable name can not contain more than 7 characters. 
(4) No special characters (such as, t+, -, etc.) are permitted. 
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Examples of FIXED POINT VARIABLE names are: I, NEMO, M5, and LEFTONE. 


How many of the following FIXED POINT VARIABLE names are not permitted in Fortran? 


INCREASE MAP KONSTANT 
N3+d 1234567T 430126 
15 JAX LEFT3 


If your answer is 3 or 4, re-read 5, and then repeat 6. 
If your answer is 5, GO TO 7. 
If your answer is 0, 1, 2, 6, 7, 8, or 9, GO TAKE A SHOWER and repeat 5. 


Examples of FLOATING POINT VARIABLE names are: 
A32 CREDIT DEBIT Y OK GAMMA FORTRAN 


How many of the following FLOATING POINT VARIABLE names are not permitted? 


NO32.6 PI ALPHjBETA 
XYZ E ALP-BET 
3.2 ALPHA 272 


If your answer is not 4, review 5 then repeat 7. 
If your answer is 4, go to 8. 


The word “FORTRAN” is a contraction of “Formula Translation". A Fortran 
arithmetic statement closely resembles a conventional mathematical formula. An 
important conceptual exception is the equal, =, sign. In Fortran, “=" means “is 
to be replaced by” rather than "is equivalent to”. Thus, in Fortran, the statement: 


Y = A3 + BIG 


means that the sum of the floating point variables, A} and BIG, is to replace the 
present value of the floating point variable, Y. 


‘Another important exception is that the left hand member of a Fortran statement 
must be a SINGLE VARIABLE NAME---it can not be an expression. 


Thus, I2 = 3 + J5 is permitted, but I2 + J5= 3 4s not. 
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The sum of 20 terms of the arithmetic series, 3 + 5 +7 ----- hi, is 440. 
This is easy, if one remembers the rule: 


Sum Total Number of Terms (Fas + Last 


2 Tern Term | 


Let us turn this rule into a formula by use of the following letters: 


s = Sa +8) 


In Fortran, the asterisk, *, is used as the mltiplication gmbol and the slash, /, 
as the division symbol. A double astrisk, **, indicates exponentiation. Otherwise, 
other symbols including parentheses are conventional. With this short briefing, 
can you write a Fortran statement for the above formula? Try it---on a separate 
sheet of paper, then 
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Congratulations if you have one of the following! Probably the second or third 
best fits the "left to right" operational philosophy. 


(A + B) * TN/2, 
(TN/2.) * (A + B) 
-5 * IN * (A + B) 
(TN * (A + B))/2. 
TN * (A + B)/2. 


nnnn om 
fl 


Note the importance of the decimal point after the 2 and before 5. Constants mst 
be floating point if the other variables of the right hand member are floating point. 
(This will be explained later.) 


FORTRAN ORDER OF OPERATIONS PHILOSOPHY 


Quantities within parentheses are computed first and the resultant value is then 
considered in the simplified expression. All operations are performed from left 
to right in the following order: 


(1) Exponentiation (**) 
(2) Multiplication (*) and Division (/) 
(3) Addition (+) and Subtration (-) 

Go To 11 
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Let us follow through an example using the correct FORTRAN "order of operations" 
philosophy. 


Example: A=B+C/D ** E * (F + 2.0) 


The calculation would proceed as follows: 


Step 1 Add F to 2.0 

Step 2 Raise D to power E 

Step 3 Divide C by the result of step 2 
Step 4 Multiply the result of step 3 by the result of step 1. 
Step 5 Add B to the result of step & 

Step 6 Replace A with the result of step 5 
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A student programmer recently wrote the following statements for the example 
below. Determine which of the statements are correct. 
Example: Add 3 to a variable, J, and square the sum. Multiply the product by 5 
and divide by 12. Let the final result equal K. 
STATEMENTS 


(4) K=5* (J + 3) ** 2/12 
(2) K=5 *J5 +3 ** 2/12 

(3) 5 * (J +3) ** 2/12 = K 
(4) K = (J +3) * 5 ** 2/12 
(5) K = (J + 3) ** 2 * 5/12 


If your answer is 1 and/or 5,’ go to 14 
If your answer is 2, 3, or 4, go to 13 


The first Law of Revision: 
Information necessitating a change in the program will be conveyed 


to the programmer after.-- and only after -- the program has been 
debugged! 


13 


2 
Statement (2), K= 5 * J + 3 ** 2/12, is incorrect. It is equivalent to: K = SJ +2, 


Statement (3), 5 * (J + 3) * 2/12 = K, is correct mathematically. However, in Fortran 
the side to the left of the equal sign can only be a single variable. Thus if the 
sides were reversed, this answer would be correct. 


Statement (4), K = (J + 3) * 5 ** 2/12, is incorrect. It is equivalent to? 


- (J +3) 5°. 
K 12 


Why not return to 12 and try that example once more. 
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The two.examples just completed involved the following statements. 
(1) A=B+C/D ** E * (F + 2.0) (2) K=5* (J + 3) ** 2/12 


In the first statement both members are in "floating mode", since they contain 
only FLOATING POINT VARIABLES and FLOATING POINT CONSTANTS. ‘In the second, both are 
in "fixed mode", since they contain only FIXED POINT VARIABLES and FIXED POINT CONSTANTS, 


The next question becomes very important----CAN THE LEFT SIDE BE OF A 
DIFFERENT MODE THAN THE RIGHT SIDE? The answer is, Yes, but only under 
the following restrictions. 


(1) The left side (FIXED or FLOATING) mst be a single variable. 

(2) If the modes of the sides differ, calculations are done using the mode 
of the right side and then put into the mode of the left. Thus in F =I + 3, 
the FIXED POINT,I + 3,is FLOATED and replaces F, However, in I = 1.5 + A, 
when the FLOATING POINT sum,i.5 + A,is FIXED, the fractional part of the 
sum is lost. Thus, if I = 1.5 + 2.23 I is replaced by 3. 
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NOW, CAN ONE MIX MODES WITHIN THE RIGHT SIDE OF A STATEMENT? 
Yes, but again only under certain cases. There are 3 such cases: 
(1) EXPONENTS 
(2) SUBSCRIPTS 
(3) ARGUMENTS OF SUBROUTINES 





The last two are described later. The first is easily learned. It seems 
logical to be permitted to use an integral FIXED POINT exponent with a FLOATING 
base. In fact, a negative base requires a FIXED POINT exponent. 

Examples: 
B = (.X)**3 
Y = TP + 3.2 - H** 


The integer, 7, is an exponent and is permitted in FIXED mode even though the right 
side is FLOATING. However, one is not restricted from using FLOATING POINT exponents 
if desired. For example: 
: Y = TP + 3.2 = H**7, 
Y = TP + 3.2 - H**3.2 
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Several Goodies have now been presented’ to you. See how many you sampled, by 


writing the answers to the following? 


(1) 
(2) 


(3) 


(4) 


(5) 
(6) 
(7) 
(8) 


(9) 


Mathematically 716. and 716 are equivalent. In Fortran, is there any difference? 
What constant does .0126 E + 4 represent? 

Three variables are IKE, MIKE, and TOM. Which one represents fractional quantities? 
Complete the limerick: 


Yes, dear Variable, the name is "Pettycur® 
It fits this data like “good” fits heaven 
But, alas I can’t use the moniker 

Since the characters number ° 


White can not equal black in color photography. Can WHITE = BLACK in Fortran? 
List the arithmetic symbols used in Fortran. 

Given? 7X + 12 = Y, write a Fortran formula which would solve this equation for X. 
Write the following as a Fortran statement and find the value of I. 


T= 2 (37 +3)— 24441)" (2 +1) ,, 
b) 


Is the statement J= 3.2 + 5.4 - 6.123 permitted in Fortran? ‘If so, what will be 
the final value of J? 
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Dear Godfrey, here are the TIPSTERS you need. I hope your’s are correct. ---~« 


(1) 
(2) 
(3) 
(4) 
(5) 
(6) 
(7) 


(8) 
(9) 


There is a difference. "716" is a FIXED POINT CONSTANT, but "716." is a FLOATING 
POINT CONSTANT, 


.0126E + 4 is another way of writing the constant 126, 


IKE and MIKE represent integers--variable names starting with I, J, K, L, M, and 
N represent FIXED POINT quantities. TOM represents a fractional quantity--it is 
a FLOATING POINT variable, since it does not begin with letters I, J, K, L, M, or Ne 


"Since the characters number more than seven! 


Yes. In Fortran WHITE = BLACK means that the value of variable, WHITE, is 

replaced by the value of variable, BLACK, 

+ add * mitiply = equal 

- subtract | divide ( ) grouping 
*® exponentiation 


X= (Y - 12.)/7, 
I = (37 + 3) * 2 = (2e(4 + 1)*929(1 + 1))/5 +1 and I= 6! 


Yes. The right member will be calculated in FLOATING POINT (2.477), and then 
FIXED variable, J, will be set equal to 2 (FIXED value of 2.477). 
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One often wants to represent a whole ARRAY of variables by a single variable 
name with individual terms of the array being designated by using subscripts with 
the array name. Provision is made in Fortran for one, two, or three dimensional 
arrays. The array name is followed by the subscripts separated by commas, and enclosed 
by parentheses. Thus Ay is defined as A(1)3 As,2 is defined as A(1, 2)$ and Aa,ay3 
as A(1,2,3). 


. The subscript must be either a FIXED POINT. constant, FIXED POINT variable, or 
FIXED POINT expression. For example; A(1), A(J), A(J + 3). 


Note that a FIXED POINT SUSSCRIPT can be used in a FLOATING ARRAY. For example, 
BETA(4).. Here is the second permissable use of mixed modes within the right side of 
an arithmetic statement. Thus, a FIXED POINT subscript may appear with a FLOATING 
POINT variable name. | 


Note that whereas the subscript is always FIXED, the array name itself can be FIXED 


or FLOATING. Thus I(1) is permitted as well as A(1)$ JET(2, 3) as well as BET(2, 3); 
etc. 
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The Compiler must know how many locations to reserve for arrays. This 
information is given to the Compiler by a DIMENSION statement which must precede 
the first reference to the arrays. One DIMENSION statement can indicate the 
range of more than one array. For example: 


DIMENSION A(200), IOTA(50, 50), PI(3, 3, 3») 


The above DIMENSION statement provides for a one dimensional array, A, of 
200 terms} a two dimensional array, IOTA, of 50 by 50 (2500) terms; and a three 
dimensional array, PI, of 3 by 3 by 3 (27) terms. 


Do you see what is required in a DIMENSION statement? The word, DIMENSION, 
4s followed by the array name, followed by parentheses enclosing the dimensional 
maximums for the subscript ,which are separated by commas. Commas also separate 


one array from another in the statement list. 
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Alphonse used a 20 by 30 matrix (array) to represent a system of equations. 
He wanted to add to the term located in the 5th row and the 23rd column. Alphonse 
wrote the following statements to indicate the range of the array and to increase 
the specified term of the array. | 


DIMENSION ALPHA (20, 30) 
ALPHA (5, 23) = ALPHA (5, 23) + 3. 


Did Alphonse GOOF? 


If you think he did, go to 21 
If you agree with Alphonse, go to 22 
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No, Alphonse did not goof---check the problem again, and then study the 


following example---for more practice! 


Alphonse also used an array, DATA, of 2000 terms. Various terms of this 
array were processed, depending upon the value of N. Thus if N = 3, then DATAs 
was processed; if N = 7, DATA7, was processed, etc. 


Let us assume that during part of the program, Alphonse desires to multiply 
the third term following the processed term, by 5 and put the result at X, Thus 
if N = 3, he wants to put 5(DATAs) at X. He writes the following statements: 


DIMENSION DATA (2000) 
X = DATA (N + 3) * 5. 


“Note the decimal point after the constant, 5, since "DATA" indicates 
FLOATING POINT. 
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Alphonse later set up the following statements for another problem. 
DIMENSION T(9), K(8) 


i) = 3.0 

T( I + J)= 4.2 

TI +J+7) =T(K1 + J 1) + 1,2 
K(7) = T(I + J) 


“tom 


teu 


The above statements would result in specific actions, On a separate piece 


of paper mark the following as TRUE or FALSE, 


The DIMENSION statement indicates a 9 term array (T) and an 8 term array (K). 
The value of Ts: is replaced by FLOATING POINT 3 

The value of Tz is replaced by a FLOATING POINT 4,2 

The value of Ts» is replaced by a FLOATING POINT 4.2 

The value of K7 is replaced by a FLOATING POINT 4,2 


If you think all 5 were TRUE, go to 23 
If you think at least 4 were TRUE, go to 24 
If you think at least 3 were TRUE, go to 25 
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The 5th one was FALSE! The value of Kz must be replaced by Ts as claimed, 
and Tg = 4.2 is in FLOATING POINT. But Kz must be FIXED POINT since K is one 
of the six FIXED POINT letter indicators (I, J, K, L, M, N). Therefore, K7 = Ta 
is a statement in which the mode of the left side of the statement differs from 
the mode of the right side. The mode of the left side determines the final mode, 
so in this case, the final mode must be FIXED POINT, since "K(7) is the left 
variable. The correct answer, then, is: 


K, =4& 


The fractional part of the result has been lost! 
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You are on the beam, and to show you how proud we are of your good work, we 
present the following story: 

Years ago in India, an old tea taster, Sanskrit, always confounded the experts 
by unfailingly being able to differentiate between BOMBAY, CALCUTTA, and KILMNJ teas. 
In those days it was common for unscrupulous tea makers to try to push KILMNJ tea 
on unsuspecting merchants. But old Sanskrit always could tell the difference, 

In spite of all efforts to discover his secret, Sanskrit carried it to his 
grave. On his death bed he is supposed to have muttered something like the following, 
“They FLOAT, but KILMNJ is FIXED", 

Years later, a far removed American descendent ‘of Sanskrit remembered these 
words. This distant cousin in writing a Fortran Compiler, and faced with. the 
problem of differentiating between FIXED and FLOATING point variables, noticed 
that "KILMNJ":when re-arranged had the alphabetic order of I-J-K-L-M-N, 

To this day, in Fortran, FIXED POINT variables begin with I, J, K, L, M, or N. 
Otherwise, they FLOAT. 


25 
What happened? Do you feel well? Let’s go over the statements. 


DIMENSION T(9), K(8) Reserves 9 locations for FLOATING POINT array, T, 
and 8 locations for FIXED POINT array, K. 
I=i 1 goes to I 7 
J=0 0 goes to J 
T(I + J) = 3.0 T(I + J) = T(1 + 0) = Ty is replaced by 3.0 
J=Jt+i J+1=0+12 1, goes to J 
T(I + J) = 4.2 T(I + J) = T(1 + 1) = Tz 418 replaced by 4.2 
(I+J+7) = Tl + J +7) = T(1 +1 +7) = Ts is replaced by 
(IeJ-1) + 1.2 TI + J -1) +1.2 = T(1 +1-1) 1.2 = Ty + 1.2 = 
3.0 + 1.2 = 4,2 
K(7) = T(I + J) T(I + J) = T(i + 1) = T2 = 4.2 is to go to Kr, 


But Kz is FIXED. Thus the fractional part is 
lost, and 4 goes to Kz, 


Let's start again, 
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So far you have been doing very well but you can’t afford to become over- 


confident. Let’s stop for a couple of minutes and try another set of exercises 
before moving on. Write the answers to the following: 


A. 


C, 


Write one statement which will reserve enough space for Ws, Wz, --- Wsso;, and 
also for a 50 by 50 matrix called MONO, > 

T 
A formla for CHI SQUARE is? cur’= TAMA) (D) = (B) (Cc) + 2 
Write this in Fortran. . (A +B) (A + D) (B + C) (C + D) 
Given; TETRA = P3 ** 2 + 5.2 + V5 Since everything in the right member is 
FLOATING POINT, the constant 2 must be 2.0 or 2.. Defend or attack this 
conclusion! 


May a subscript be a FIXED POINT constant? If yes, give an example. 

May a subscript be a FIXED POINT variable? If yes, give an example. 

May a subscript be a FIXED POINT expression? If yes, give an example. 

May a subscript be a FLOATING POINT constant, variable, or expression? 

Maya subscripted variable name be FIXED or FLOATING? If yes, give an example. 
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2? 
CHECK YOUR ANSWERS! 
DIMENSION W(350), MONO (50, 50) 
CHI=T* (A*D-B*C+T/2,) * 2/((A +B) * (A +D) * (B+ C) * (C + D)) 


It is not necessary to give the decimal point! Nor is it wrong ! Keep in mind, 
that an exponent may be a FIXED POINT quantity even though the base may be 


“FLOATING POINT, This is one of the cases where a FIXED quantity is allowed 


to be mixed with FLOATING variables. A second case is SUBSCRIPTS, The third 
case, which you will get later, involves the uses of ARGUMENTS of functions. 


Yes, a subscript can be a FIXED POINT constant; ex. U(3) 

Yes, a subscript can be a FIXED POINT variable; ex, U(NEND) 

Yes, a subscript can be a FIXED POINT expression; ex. U(I + J - 3) 

No, a subscript can not be a FLOATING POINT constant, variable, or expression. 
Yes, the subscripted variable name can be FIXED or FLOATING, ex. J(3) or T(3) 


If you had significant differences from the above, REVIEW 18 and 19, then 
Go To 28 





A Fortran 80 Column Statement Card 
Note Columns 1 - 6, 7 = 72, 73 - 80 





5 
CONTINUATION INDICATION col.6 
FORTRAN STATEMENT cols. 7 « 72 
SERIAL IDENTIFIERS cols. 73 - 80 
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You just saw a FORTRAN STATEMENT CARD. Did you really GRASP the important 
facts presented? Wait! Don’t look at it! Try these questions first! 

In what column did the Fortran statement start? 

In what column must the Fortran statement end? 

In what columns are the Statement Numbers given? 

If there is a CONTINUATION indication---what column should it be in? 

Where are Serial Identifications made? 


Now, look at 28 again! 


If you feel you did 0.K., then 
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Column 7 is the reference column to remember, when writing Fortran programs! All 
statements or continued statements begin there! Starting at column 7 statements 
progress through column 72 --- ONE CHARACTER PER OOLUMN, Remember, only ONE statement 
is allowed on the same card. 


If a statement extends beyond column 72, then a CONTINUATION card is required. Column 6 
is used when CONTINUATION cards become necessary. A blank in column 6 indicates no 
CONTINUATION; NINE digits 1, 2, 3, --- 9 indicate corresponding CONTINUATIONS., Thus, 

if one statement requires 3 cards; the first card has a blank in column 6; the second 
card has a 1 in column 6; and the third card, a 2 in column 6, Remember, each continued 
statement also starts in column 7 -- following the continuation digit in colum 6, 


Statement numbers are given in columns 1 - 5, All statements need not have a number-- 
only those which are later referenced in the program by other statements. 


Columns 73 - 80 need not be used as the compiler does not process these. However, 
they may be used for serializing the cards, or providing other information. 
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COMMENT CARD 
Note "C" in column 1 
Comments start in col. 7 






The "C" in col. 1 tells 
the compiler that this is not a 
Fortran statement card and the data 
starting in column 7 is not processed. 
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You have seen the 2 types of cards which are source programs prepared for input 
to the compiler--STATEMENT cards and COMMENTS cards. A key punch operator will 
probably punch these cards for you. However, it is important for you to know card 
format so that you can write programs in correct form for the operator to punch-- 
otherwise, they are liable to punch you! Your programming will probably be done 
on printed forms similar to 33. (These FORTRAN CODING forms are available from 
the Publications Dept.) 


33 portrays the requirements of card format you have just reviewed. Study 
this example carefully, and note the column sub-divisions. 


P. S. We are assuming that you have “picked up" the correct concept 
that all constants used in Fortran--including Statement Numbers-- 
are our "good old" decimal numbers--the Compiler takes care of 


all conversion’s for you. 


Go To 33 





1604 FORTRAN CODING FORM NAME. 
PROGRAM Doomed PAGE 1 
ROUTINE Shot to the Moon DATE 6/7/99 






FORTRAN STATEMENT 










PROGRAM DOOMED 







Cc A PROGRAM TO END ALL PROGRAMS 125 
DEATH = SHOT3 + 5.2 | 126 
CHI = T * (A*B ~ C*D + T/2.)**2/( 127 
(A + B) * (C + D).* (A + D) * (B + C)) 128 
GO TO 6 129 







Serial 124 is a PROGRAM NAME CARD -- required for identification 
Serial 125 is a COMMENT CARD (Indication in Column 1) 

Serial 128 is a CONTINUATION CARD (Indication in Column 6) 
Serial 129 is a CONTROL STATEMENT CARD (To be discussed later) 


Go To 34 


a 


How well do you know Card Format? Try matching the following numbers and 
letters on a work sheet. 


1. Serialized Card numbers A. True 
2. Column 6 B. Punch in columns 1 and 6 
3. Norml Statement Card C. Col. 1 is blank or numeric 
&u "C® in column 1 D. Yes 
5. Card spaces are ignored by the compiler E, Two Fortran Statements 
6. Can a. statement of more than F, False 
660 characters be used? G. Comments Card 
7- Gan not be put on one card. H. Continuation indication column 
S.A Fortran statenant shoul 1. tua 73 
« No 


9. Statement numbers need not be 
given in numerical order. 
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The correct matching should be: 
(1I, 2H, 3C, 4G, 5A or SD, 6J, 7E, 8A or &D, 9A or 9D) 


If you missed a couple, blame us since the test contained a few "sneaky items", 


Are spaces ignored by the compiler? Yes. We hoped you would assume this to be 
true. Note, in example of 33, we used some spaces. 


Is 660 the maximum number of characters permitted? Yes. There may be 9 
CONTINUATION CARDS plus the original card--making 10 cards total. Each card uses 
columns 7 through 72 (66 characters) for a grand total of 10(66) = 660. Note, the 
660 means 660 columns including any spaces that are involved within the statement. 


Also statement numbers need not be in numerical order--we had not mentioned 
this previously---but now you know it! 


Go To 3% 
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The name of a program is "ANIMAL", Ten variables: A, ALPHA, B, C, D, HIPPO, 
LYNX, BEAR, I, and J are involved. 


HIPPO, LYNX, and BEAR are subscripted. HIPPO is a 60 by 60 matrix, LYNX is an 
array of 12 terms, and BEAR is a 16 ‘by 16 by 16 matrix. The program is identified 
as the "ANIMAL PROGRAM". The following formula is used in the program and is 
referred to later as statement 22. 


HIPPO,, so = LYNXg, + (I + J)? 


Also the following formula is used and is subdivided into CONTINUATION statements. 
For this example, assume the continuations start with each slash (division) shown. 


BEAR, = A(C + D)* / B + (B) (C + D)°/ (ALPHA + 3.2) 


Program the above example, using correct card format, and indicating the card 
column sub-divisions. 


Go To 37 
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FORTRAN STATEMENT 
COLUMNS 















PROGRAM ANIMAL 
ANIMAL PROGRAM 

DIMENSION HIPPO (60, 60), LYNX(12), BEAR(16, 16, 16) 
HIPPO(1,50) = LYNX(11) + (I # J) ** 3 

BEAR (I) =A * (C+D) ** 2 

{B+B* (C +D)'** 3 

/ (ALPHA + 3.2) 


Go To 38 
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A computer program is similar to the recipe for "making a cake” in that certain 
cake ingredients are listed followed by a set of directions or order of events to 
put these together. - 


Up to this point, we hope you have learned the Fortran ingredients: symbols, 
card format, rules of constants, variables, etc. We now need to learn certain 
CONTROL statements which will give direction to the normal sequential order of 
events of a written program. 


Although Fortran includes a few more CONTROL statements than will be presented 
in the course, the following seven will be sufficient to provide efficient control 
for any Fortran progran. 


GO TO (unconditional) CONTINUE 
GO TO (computed) PAUSE 
IF 

DO STOP 


Go To 39 
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The first CONTROL statement is the unconditional "GO TO". This statement has 
already been used at the bottom of each page -- where the number expressed is a page 
number instead of a statement number. Thus, the correct form is, GO TO n, where n 
is a statement number. An example is: GO TO 372. . 


The second CONTROL statement is similar. The computed "GO TO” must include a | 
list of statement numbers, separated by commas and enclosed by parentheses. This 


list is followed by a comma and a non-subscripted fixed point variable, such as I. 
The value of I then determines which branch to take. 


For example GOTO (100, 72, 136), J 
if J =1, control will go to statement 100 
if J = 2, control will go to statement 72 
if J = 3, - control will go to statement 136 
if J = 0, 4, or larger (in this example), an error halts 
the program. 


Go To 40 


My) 


Although the first GO TO is self evident, the second needs some practice. 
Below are some further examples of computed GO TO statements. 


GO TO (45, 129, 5, 5, 5),JAKE 

If JAKE = 1, GO TO Statement Number 45 

If JAKE = 2, GO TO Statement Number 129 

If JAKE = 3, 4 or 5, GO TO Statement Number 5 


GO TO (7, 18),M 
IF Mis 1, GO TO Statement Number 7 
IF M is 2, GO TO Statement Number 18 


You are to be guided by the GO TO Statement below, assuming NN is 3. 


Go To (40, 40, 41)NN 


R, U, HESITANT: born 1782, before advent of modern computer; known 
as "father of the flow chart", planned events of each day with a 
diagram similar to following: 







Yes End of Day? 
(1) See Indian? 


Look 
No Again @ 


44 
Given the following flow chart of part of a program, write Fortran statements 
which will do everything shown in flow chart. 








SUMA is a subscripted 
variable (9 by 9) 
DATA is an array of 81 
terms | 











Programmer comments 
“SYSTEMS ERROR SUM*® 





Go To 42 
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Your Fortran program from the flow chart should be the following* 


Columns 
i ---5 6 7 
PROGRAM ERROR 
Cc SYSTEMS ERROR SUM 
DIMENSION SUMA (9, 9), DATA (81) 
I=1 
J=1 
3 SUMA (I, J) = DATA (J) ** 3 + 3.2 
h IT=I+1 
J=J+1 
GO TO (3, 3, 3, 5)si 
5 SUMA (I, J) = DATA (J) ** 2 + 3.2 
GO TO 40 


If you missed some, REVIEW 38-40 and try 41 again 
Otherwise: 


Go To 43 
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The next CONTROL statement is the "IF" statement.. 
IF (X - Y) 20, 30, 15 


This example means, "If the value of X - Y is less than zero, go to statement 20; 
if equal to zero, go to statement 30; and if greater than zero, go to statement 15." 


IF (X) 20, 30, 15 


This example is identical except the value of "X" 4s now examined as to whether 
it is less than, equal to, or more than zero. 


IF (J + 2) 15, 15, 20 


In this example, "J + 2" is examined. If J + 2 is less than or equal to zero, 
control goes to statement 15; if more than zero, to statement 20. 


Go To 44 


by 


Summarizing the IF statement then the format is: 


FIXED or FLOATING VARIABLE Less than Equal to Greater than 
or _ zero branch zero branch zero branch 
FIXED or FLOATING EXPRESSION Statement Statement Statement 
number , number, number 


Following the parentheses are 3 statement numbers separated by commas. These 
are the paths or branches. If the quantity in parentheses is less than zero 
control goes to ‘the first statement number, if equal to zero, control goes to the 
second statement number, and if greater than zero, control goes to the third 
statement number. All 3 statement numbers must appear. 


Another example; 
IF (BETA(I)) 10, 20, 10 


Go To 45 


45 
Write IF statements for the following three examples. 
(1) The variable JET is usually negative, However, it can become zero. If JET 
ever becomes greater than zero, an error results. As long as JET is negative, 
the branch is to statement 115; if zero, it is to statement 22. The error stop. 


is at statement 5. 


(2) If the expression "DATA + Y - .5E2" is zero, the branch is to statement number 6; 
otherwise, the branch is to statement number 100. 


(3) If W* is negative the control goes to an error exit at statement 42: if zero, 
control goes to statement 5, otherwise to statement 15. 


Go To 46 
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Answers to the 3 "IF" examples are: 
(1) IF (JET) 11, 22, 5 

(2) IF (DATA +Y- .5E2) 100, 6, 100 
(3) IF (W** 2)" ha, 5, 15 


Did that constant in example 2 look strange? Remember, .5E2 is a FLOATING 
POINT constant equivalent to .5(10)? = 50. 


Note in second example, the IF statement becomes a two way branch. Also note 


the "pseudo example” given below, is a_one way branch or an unconditional GO TO, 





IF (YOU MISSED ANY OF ABOVE) 43, 43, 43 


Otherwise, 


Go To 47 
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"How do I get to Casey street?", I asked the cop. 

"Go to 15th, and check the SIGNAL there. If red, go right to 50th, if yellow 
turn left to 34th, and if green, go straight ahead to 11th." 

"But I’m color blind", I protested, “how can I tell if the signal is red, 
yellow or green?" 

“It’s easy", said the cop. “When you get there check the amount of money you 
have in your left coat pocket. If none, the light is yellow, if at least one cent, 
the light is green, if you find a hole in your pocket, the light is red." "Very 
simple", he concluded boastfully, “if you know FORTRAN", 

I thought for a minute, wrote the following on a sheet of paper, am handed it to 
the cop, with the question, "Is this what you mean"? 

GO TO 15 
15 IF (MONEY) 50 34 11 

"That’s right” said the cop, "except you forgot the comms between the statement 

numbers." 


Go To 48 
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Write a program for the following flow chart: 


A and DELTA 
are each & by 5 
arrays; BETA is 
a 4& term array 


Set DELTA 
i,j 
to ZERO 


3 ‘BETA, / 











DELTA, |; A 


LESS THAN (10) 


EQUAL 
REATER————>(50) 


Go To 49 





LESS THAN, EQUAL TO, 
QR GREATER THAN ZERO 





Reprinted from the abstract of SECOND ANNUAL CONFERENCE ON PROGRAMMING 
ETHICS, Washington, D. C., 1961. 


The time has come for programmers to: 


STOP adopting programs, written by other people, for unrelated problems, 
STOP using facts when one can use scientific truthse 

STOP demanding computer operations with extra sensory perception abilities, 
STOP blaming the computer for writing bugs into the programe 

STOP trying to make complicated problems out of simple ones. 


In other words ----- STOP PROGRAMMING! 


10 
11 


12 
20 


bg 


Fortran program for the preceding example should resemble the following: 


PROGRAM TEST 

TEST RUN NUMBER 100 

DIMENSION A(4, 5), DELTA(4, 5), BETA(4) 
I=1 

J=1 

IF (A(I, J)) 11, 12, 11 

DELTA(I, J) = BETA(I)/A(I, J) 

GO TO 20 

DELTA(I, J) = 0. 

IF (DELTA(I, J) - 6.2) 10, 30, 50 


If your program differed significantly, review 43 through 47, and do 48 again. 


Go To 50 
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Most Fortran programs involve looping, the repetition of groups of statements 
with varying parameters. The control statement involved is the "DO" statement. We 
begin with a very simple example. 


K=0 
DO 303 J=1, 3 
K=K+2 


303. SUM(J) = J ** 2 +K ** 2 


The above DO statement calls for 3 loops, since the loop control sequencer 
®J= 1,3" indicates an initial value"l", and a terminal value,"3." The statements 
involved in the loop are all the statements following the DO down to and including 
the statement number 303. Thus: 


First loop K=0+2 

SUM, = 1° + 27 
Second loop K=2+2 

SUM2 = 2° + 4? 
Third loop K=4%+2 


SUM, = 3° + 6? 
Go To 51 


How about another example? 


DO 7 NN=1, 4 
K(NN) = 3 ** NN + 5 
L(NN) = NN ** & +5 

7 (NN) = K(NN) + L(NN) 


The loop control sequencer "NN = 1, 4" indicates 4 loops since the initial 
and terminal values are 1 and 4 respectively. The "7" indicates that all statements 
following the DO statement, down to and including statement 7 are indicated in the 
loop. 


Ky = 3° + 5 Ky = 3°? + 5 
Loop 1 0 la = 1" + 5 Loop 3+ Ls = 3° + 5 
M, = Ka + La = 14 Ms = Ks + Ls = 118 
Ke = 32 +5 Ky = 3° +5 
Loop 2 La = 2° +5 Loop 4 Le = 4° +5 
M2 = Kz + La = 35 My = Ky + Le = 347 


Go To 52 


A. 


B. 
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Write the DO statement which will include the following statements in a loop 
using an initial value I = 1 and terminal value I = 200. 


BETA(I) = ALPHA(I) -D** 2+7.0 
IF(BETA (I)) 10, 11, 16 
16 ALPHA(I) = F(I) - 4.2 


Given the DO statement below, answer the following questions. 
DO 100 JK = 5, 27 
1. 


The "100" indicates that the loop should be executed 100 times (TRUE, FALSE) 
2. 


The control sequencer, JK = 5, 27, indicates which one of the following: 
(1) JK takes on values 5, 6, 7, ---- 27. 


(2) J takes on values 1, 2, ---- 5 ami K takes values 1, 2, --- 27. 
(3) J becomes $ and K becomes 27. 


Go To 53 
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For the first example A, of the preceding card, you should have the statement-- 
DO 16 I = 1,200 


For the second example B, the correct answers are; 
1. FALSE 


2. JK takes values 5, 6, 7 «== 27 


If you had some trouble with these two examples, GO TO 50. Otherwise, 


Go To 54 


5+ 


Lest you become over-confident, let us warn you that there is more than meets 
the eye in the DO statement--that is more than has met your eye to this point! 


The first of these “eye hiders" is the fact that the control sequencer may 
contain a third variable or constant, following the initial and terminal values. 
For example, 





DO 17 I=1, 51, 5 


Here, the control sequencer, I = 1, 51, 5, indicates the following: 


1 is to be the INITIAL value of I 
51 is to be the TERMINAL value of I 
5 is to be the INCREMENTAL value of I 


When the INCREMENT is equal to 1 (the usual case) this third variable or constant 


need not be expressed. Thus I takes values 1, 6, 11, ------ 51, since the 
incremental value in the above example was 5. 





Go To 55 
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The index of the "control sequencer® must be a fixed point variable. That is, 
the index must start with one of the letters I, J, K, L, M, N. (Actually, the 
compiled program works more efficiently if the single letters I, J, K, L, M, N are 
used for the index). 





Of utmost importance is the fact that the two (or three) parameters (the initial 
value, the terminal value, and the increment) following the index of the control 
sequencer can be FIXED POINT VARIABLES as well as FIXED POINT CONSTANTS but not 
EXPRESSIONS. Side by side examples are: 


DO 22 I=1, 10 DO 22 I=N, 10 
DO 56 J=4, 100, 2 DO 56 J=N, K, 2 
DO 30 L= 5, 89, 6 DO 30 L=N, K, J 


Go To 56 
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There is still more to learn about the DO statement! For example, 


The loop itself involves learning a few concepts. 

One of these has to do with TERMINAL values of the index. The index is 
increased by the INCREMENT value at the bottom of the loop--after the last statement 
of the loop. If the index variable becomes greater than the terminal value after 
being incremented, the DO loop is terminated, and the index value is set to zero. 
The program then goes to the next sequential statement following the last statement 
of the DO loop. If a DO loop is left before the loop is terminated, the last value 
of the incremented index is retained. 


For example: 
D017 J=3, 7, 2 


LOOPING continues, for three loops (J = 3, J = 5, J = 7) and on normal exit 
J is set back to zero. However, if a inner exit (programmed) had occurred during 
the second loop, the last value, 5, would have been retained by J. 


Go To 57 
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Diagramatically, we might show the looping concept as follows: 


Initial Index Value 
Terminal Statement Number 










Terminal Index Value 
Increment 


Index 





DdO71 I=1, 7, 3 
ALP(I) = BBTA(I) + 8.9 
J=J+i1 

71 ‘BETA(I) = BETA(I +1) + J 


'rsi+3 | 


Next Statement in Program 


Go To 58 


58 


There are two very important restrictions concerning the last statement of a 
DO loop. First, there must be a statement number (in columns 1 ~ 5) for this 
statement. Second, the last statement can not be a control transfer type statement 
("GO TO" or "IF"), 





However, there are instances where you will want to put an IF or GO TO 
statement at the end of your DO loop. If so, a special "do nothing” type statement, 
CONTINUE, is required to end the loop. Since this is the last statement of the 
loop, it must have the same statement number previously expressed in the DO 
statement. An example follows: 


DO 6 J = 5, 100 
B(J) = A(J) + TIME ** 3 
IF (B(J)) 10, 6, 6 

6 CONTINUE 


Go To 59 


59 


Finally, DO loops can be contained within DO loops as long as the inner DO 
Loops are nested loops (contained within the outer loop). It is possible under 
this condition for several DO loops to terminate with the same statement. Also 
note that nested DO statements can not use identical index variables. The following 
diagrams indicate correct inner loops. 











DO 31 =1, 100 DO 31 =1, 100 
| DO 4K =1, 100 DO 3K = 1, 100 
Ty 

DO 5J=1, 100 DO 3 J = 1, 100 


Go To 60 
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Now that all the rules are in, did we over DO it? The best way to find out is 
by examples such as the following. 


A programmer recently wrote the following portion of a problen. 


PROGRAM TEMP 

DIMENSION W(500), 2(100), F(100) 
DO 15 I = 1, 100 
F(I) = DATA + 7.0 

12 G=G+t. 
W(I) = F(I) + HL3 ** 2 - (Z(I) - G) 
IF (W(I)) 12, 12, 15 

15 CONTINUE 


Did this programmer violate a DO loop principle? If you think he did, Go To 62. 
If you think the program is correct, Go To 61. 
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The statements were OK. It is good that you noticed that a CONTINUE statement 


must follow a transfer control statement when the transfer 


The following incident took place one day: 













“Hey Moe, I just 
wrote this DO loop 
and it looks good 


N, 100 


Gal 


L+ 3 


If you agree with Joe, go to 63 
If you agree with Moe, go to 64 


ends a DO loop. 









"No Joe, it looks 
like something is 
missing but I don’t 
know what it is." 








end of the DO loop: 
of the loop. 











YOU 

JUST TOOK 
A 

DETOUR 


Since the progra 
was correct 


Notice that it was desired to have a control transfer (IF statement) at the 
Thus it was necessary to include CONTINUE as the last statement 


Go To 61 
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Roses are red, 
Violets are blue; 
Joe was WRONG 
and so were YOU. 


The last statement of a DO loop without a statement number is like a dog 
without a tail---something is missing! 


Go To 64 


Moe was right--and so were you. 


THE STATEMENT NUMBER OF THE LAST STATEMENT OF A DO LOOP CAN NOT BE MISSING! 
There are two errors in the following example. Fird then. 


DO 32Il=J, 3 +5 
32 N(I) = M(I) -1I 

DO 14 K = 1, 10 

DO 16 L = 1, 10 
14%. GAMMA(K) = D(K) + F 
16 D(L) = F ~ 5.2 


Go To 65 
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The two errors in the preceding example were: 


1. The TERMINAL value parameter of a DO loop cannot be an expression 
(J + 5)-<-it is limited to FIXED POINT VARIABLES or FIXED POINT 
CONSTANTS. 


2. The inner DO loop (WO 16 L = 1, 10) is not nested since statement 
16 ig not included within the outer DO loop (DO 14 K = 1, 10). 


If you missed (1) above, REVIEW 55, then 


If you missed (2) above, REVIEW 59, ven——} > Go To 66 
If you had then both correct (teacher’s. pet ) 
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Some goofs are self evident--yet easy to make. Following are examples. 


DO loop index 
DO 12 A= 1, 10 must be a FIXED 
t+___—___— | POINT VARIABLE 


Index of inner DO loop can not 


be same as index of outer loop. 


ee DO 20 I= 1, 10 
20 BETA(I) = GAMMA(I) + 7.0 
4g  ALPHA(I) = BETA(I) + 50 


Go To 67 
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Three control statements remain. These are "duck soup” for "pros" who have 
progressed this far! ; 


First, there is the PAUSE statement, Example. PAUSE 712 

This will cause the program to PAUSE, indicate the PAUSE number, ard permits the 
operator to make one of the following decisions: (1) abandon the program, (2) proceed 
with no manual changes. The number following the PAUSE (up to 5 digits are permitted) 
enables the PAUSE statement to be uniquely identified. By using PAUSE statements, 

one can periodically determine where he is in the program. 


Second, there is the STOP statement. Examples STOP 355 

This causes the program to be abandoned after indicating the STOP number. The 
number (up to 5 digits) provides identification in case one has a STOP at various 
places in the program. 

Third, there is a REWIND statement. Example? REWIND 3 

This will cause magnetic tape number 3 to be rewound to its starting position. 


Go To 68 
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The next few cards will attempt to lead you through a General Review. Good luck! 
Write statements for the following Flow Chart. 


Flow Chart of “TEMPEST® Program 














Is CG, less than, 
equal to, greater 
than zero, where 

i =1, 2, ... 400 





g 
terminated? 
m6 


Go To 69 
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How does your version stack up with the following? 


PROGRAM ‘TEMPEST 

DIMENSION A(400), B(400), C(400) 

K =0 

pO 12 I = 1, 400 

C(I) = A(I) ** 3 + BCI) ** 3 

IF (C(I)) 3, 4, 12 . 

C(I) = QO. band C(I) 

GO TO 12 | 
& K=K+1 : 

' IF(K - 100) 12, 50, 12 

12 CONTINUE 

STOP 400 
50 STOP 100 


You are permitted some slight variations from the above. For example, K = 0 
can appear before the DIMENSION statement. Also, Statement 50 can be placed 
4nside the DO loop. After reviewing your errors, 





Go To 70 
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The following frequency table, containing FIXED POINT CONSTANTS, lists 100 
corresponding values of K and L. 





MEAN and VARIANCE are found using the following 
rules, where all summtions range from i = i, 


to i = 100. 
MEAN 2% — 
2 _( 
vamavce core Ea)" 


N-1 


where N = > L, 


Write the program, "MEASURE", to find MEAN and VARIANCE in FIXED POINT, 


Go To 71 
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Program “MEASURE” would appear similar to the following: 
PROGRAM MEASURE 
DIMENSION K(100), L(100) 


JVARNCE = 0 
DO 5 I = 1, 100 
N= L(I) +N 
MEAN = K(I) * L(I) + MEAN 
5  JVARNCE = L(I) * K(I) ** 2 + JVARNCE 
MEAN2 = MEAN ** 2 
MEAN = MEAN/N 
JVARNCE = (JVARNCE - MEAN2/N)/(N - 1) 


100 
Note MEAN2 is assigned to( > X,L, )*; fixed point VARIANCE = JVARNCE 
i= 1 


Study the above! You may have had another correct version. If so, and if it is 
more efficient, please send a copy to Senior Senator Humphrey, St. Paul 1, Minnesota, 
Go To 72 


A sine function is evaluated by using a truncated Chebyshev Polynomial of 
the following form: 


SINE X = Cy X + C3X® + CgX* + CrX” + CoX® + CxaX** 
where X is the argument in radians, and Ca, C3, --- Cii are Floating Point Constants. 


Assuming 200 arguments located in array, X, and the constants Cy, C3 o-- Cy 
located in the 6 member array C, write a Fortran program to compute 200 SINE values. 


Go To 73 
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Two solutions are given. The first uses a nested form of the polynomial. 


PROGRAM SINES 
C USING CHEBYSHEV SERIES NESTED 
DIMENSION X(200), C(6), SINE(200) 
DO 10 I = 1, 200 
= X(I) ** 2 
10 SINE(I) = X(I)*(C(1)+5*(C(2)+4S*(C(3)4+8*(C(4)+5*(C(5)+5*C(6)))))) 


Using the polynomial in its original form provides another solution. 


PROGRAM SINES 

C USING CHEBYSHEV SERIES NOT NESTED 
DIMENSION X(200), C(6), SINE(200) 
dO i0 T= 3 200 





SINE(I) = 
L=0 
DO 10 K = 1, 11, 2 
L=L+1 
10 SINE(I) = SINE(I) + C(L) * X(I) ** K 


As you see, theres more than one way to "skin the sine". Study the above and then 


Go To 74 
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It was a dismal day--even the local pub, down the street, looked deserted. I 
paced nervously back and forth--there was so little time to save Detson! 


"My dear Watson", said Holmes, looking yp from the clues of the Detson case, 
“quit your pacing and read this note taken from the body of Lippy Fug." 


Dear Holmes, 


Reading this means you finally got me! But you’re a gentleman, Holmes, and 
I always give a gentleman a break. Calling a certain number can save Detson-- 
and earn you his fee. So I leave a puzzler, me lad; The sum of all the possible 
products, taken two at a time, of all the convict numbers housed at Ole Bailey 
and Ole Gimey is first needed. This sum divided by 12, rings the right bell, 
old chappie. 


"You see, my dear Watson", continued Holmes, “I called the warden and found 
Ole Bailey has 100 convicts numbered 50 to 149; Ole Gimey contains 100 convicts 
numbered 1 to 100*. 


"But Holmes, I shouted, you know am I know that 100 by 100 gives 10,000 
possible products which must be summed and then divided by 12. Detson dies in two 
hours, we can’t possibly do it by then. 


"Calm yourself, dear Watson", said Holmes. “I just called the firm of Control 
Data, Ltd. They have a gadget called the 1604 which can get us the answer. However 
the problem must be written in Fortran language first. They told me a student now 


learning Fortran should be able to write this program. I am waiting for his call 
now, and as soon as he finishes, we will all 


Go To 75 
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PROGRAM DETSON 

CALL NUMBER LOCATED AT RING 
DIMENSION BALLEY(100), GIMEY(100) 
RING = 0. 

GENERATE OLE BAILEY CONVICT NUMBERS 
BAILEY(1) = 50. 

DO 20 I = 2, 100 

BAILEY(I) = BAILEY(I - 1) + 1. 
GENERATE GIMEY CONVICT NUMBERS 


GIMEY(1) = 1. 
DO 30 I= 2, 100 
GIMEY(I) = GIMEY(I - 1) + 1. 


DO 40 I = 1, 100 
DO 40 J = 1, 100 
RING = RING + BAILEY(I) * GIMEY(J) 
RING = RING/12. 
STOP 
Go To 76 


(Note? We have inserted 
two comments to indicate 
what the program is doing 
at certain times) 
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The Detson case ended not only Lippy Fug’s secret, it also ended the review 
of the first phase of this course (definitions, symbols, and control statements). 


Incidentally if you objected to the Detson case as being unrealistic because 
telephones were not available in Sherlock Holmes’ day, let.us point out that 1604's 
by Control Data Corporation were also unavailable. 


The point we want to make iss: If you had trouble with any of the review 
exercises up to this point---REVIEW! REVIEW! REVIEW! Then try the exercises again, 
If you then feel like Floyd Patterson, after his first fight with “Ingy", REVIEW 
some more. 


If you feel like Floyd after his secon) meeting with Ingy, 


Go To 77 
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As the art of programming developed, it became evident that certain mathematical 
functions (for example; sine, cosine, logarithm, square root, etc.) were used several 
times in the same program or they were used in many programs. Programmers soon realized 
that much time and inconvenience could be saved if these functions were written as 
packaged routines which could be referenced by main programs when desired. These 
packaged routines became known as SUB-ROUTINES, 


As the number of these SUB-ROUTINES increased, it was a common procedure to store 
them on some on-line peripheral input. device and to bring them into the computer when 
they were referenced by a program. Magnetic tape units were particularly effective 
in this respect. . 


The 1604 Fortran System uses such a tape and it is known as the REFERENCE 
LIBRARY. Tape Unit 1 is reserved for the REFERENCE LIBRARY. This tape contains 
mathematical SUBROUTINES which are always available for the programmer to use. 


~ 
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The following SUB-ROUTINES are available from the REFERENCE LIBRARY. 


ACOSF - - - - - arc cosine ABSF -.-.- absolute value 
ASINF - --- - arc sine EXPF - . . = - "e" to a power 
ATANF - - - - o arc tangent SINH - - - - - hyperbolic sine 
SINF - - - -- sine COSH - - - - - hyperbolic cosine 
COSF ----- cosine TANH - - - - - hyperbolic tangent 
TANF ----=- tangent LOGF ----- natural log 

SQRTF - - - - - sq. root LOGIOF = - = = common log 


These are referenced as shown in the following statement examples: 


Example 1 Y = SINF(X) + 7.7 
Name of Argument must be in 
Sub-routine parentheses; and must 
(must be one . be either a FLOATING 
of the forms POINT: VARIABLE, CONSTANT, 
given above) ‘or EXPRESSION, 
Example 2 ALPHA = BETA(I) + LOGF(X + Y) 
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The following examples and their solutions should help. 


Given? Ax* + Bx +C=0. Solve for positive x. 
= (.B + SQRTF (B**2, - 4 *A*C))/(2.*A) 


Given: Secant B + sind = Y. Write in Fortran. 
Y = 1./COSF(BETA) + SINF(ALPHA) 


Given: G = log, (x+ 4% 2-1 +), Write in Fortran. 
G = LOGF(X + SQRTF(X**2 - 1.)) 
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Subroutines are often needed which are not in the REFERENCE LIBRARY, or you may 
want to sub-divide your program into sub-programs. In either case, you will need to 


know how to write subroutines, where to locate them in the overall program, and how 
to call them up when needed. 


Let’s first learn how to write a SUBROUTINE! One simply writes an identification 
statement (see'a’ below) which names the SUBROUTINE and lists its arguments. Following 
this are the conventional statements (which you have learned) and these make up the 
body of the SUBROUTINE (lines below). Within the SUBROUTINE, you may desire a RETURN 
statement (see b). A RETURN statement will provide an inner exit from the SUBROUTINE.- 
before reaching the normal (regular) exit. The normal exit is formed by writing one 
END statement (see c), as the last statement of the SUBROUTINE. Study the example. 

















a. SUBROUTINE BETSY (ALPHA, NET) 
be RETURN 
Ce END 
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Having written a SUBROUTINE, where does one place it in the overall program? 
The correct order of placement is: PROGRAM NAME, main program containing conventional 
statements, the main program END statement, first SUBROUTINE with its END statement, 
second SUBROUTINE with its END statement, the next SUBROUTINE with its END statement, 
etc. Finally, an extra END statement, as a final statement, indicates the conclusion 
of the overall progran. 


Every calendar lists a correct order of months and days. Thus it would hardly 
be "kosher" to try to place June after October, Fortran is like this in that the 
order above must be followed. 


Lest you forget, we re-emphasize this order by having you 
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The diagram indicates the organization of a complete example program. 


PROGRAM EXAMPLE ————) Program NAME Statement 





Main Program _ 








Statements_ 
‘SUBROUTINE ONE(X, Y) 
__Subroutine _ 
__Statements _ 
Statement 
Sub-Routine 
5 TWO(G, ALPHA) Statements. 
—Subroutine 
Statements 
END LAST Statement 
END ——________} of 
Program 
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One phase of SUBROUTINE usage remains -- how to reference (CALL) the SUBROUTINE 
from the main program! This is done by a CALL statement which specifies the name and 
arguments of the SUBROUTINE. After the subroutine is executed, control goes to the 
statement that follows the CALL statement. 


Of particular importance are two requirements: first, the number of arguments 
in the CALL statement and the identification statement of the SUBROUTINE must be equal; 
second, the modes (FIXED or FLOATING) of the arguments in the CALL statement and the 
corresponding arguments of the SUBROUTINE identification statement must be the same. 


For example: 
CALL GENER (ALPHA, BETA, JET) 


SUBROUTINE GENER (X, Y, N) 
Note: 
Each contains 3 arguments and corresponding arguments agree in mode. 
The SUBROUTINE arguments (X, Y, N in this example) are DUMMY VARIABLES.-they 
always take on values of the: corresponding arguments in the CALL statement. 
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Can you write a program involving 












Find Y using 
SUBROUTINE ROOT 
with arguments 
D, E, F, YX 










Arguments 





rol 
a SUBROUTINE? Use the flow chart below: 







Is Y 
larger than 
zero? 








No GAMMA = 
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PROGRAM OUTLINE 
CALL ROOT (D, E, F, Y) 
IF (Y) 10, 10, 20 
10 GAMMA = ALPHA/SQRTF(1. - Y) 
GO TO 30 
20 GAMMA = ALPHA/SQRTF(Y) 
30 END 
SUBROUTINE ROOT (A, B, C, X) 
IF (B**2-4 *A*C) 1, 2, 3 
1 X=0. 
RETURN 
2 X = (-B)/(2. * A) 


RETURN 

3 X= ((-B) + SQRIF (B ** 2-4, *A*C)) / (2. * A) 
END 
END 


How did you do? If you had trouble, try reviewing 80-84. 
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Two arrays A and B eaoh contain 100 terms. A subroutine, TEST, examines each 
set, made up of two corresponding terms (one from each array), and determines which is 
the larger. If the term from A is larger than the corresponding term from B, the 
subroutine assigns the digit 1 to J. If the term from A is smaller, the subroutine 
assigns the digit 0 to J. If the terms are equal, B is assumed the larger. 


The main program, TESTER, always computes the following: 


X = sine of the larger term of a set 
Y = the square root of the smaller term of a set 
Z = the square root of the larger + the cube of the smaller. 


Write the program, TESTER, including the subroutine, TEST, which will determine 
X, Y and Z for each of the 100 sets. 
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PROGRAM TESTER 
DIMENSION A(100), B(100), X(100), ¥(100), 2(100) 
DO 15 I = 1, 100 
CALL TEST (A(I), B(I), K) 
IF (K) 14, 14, 12 
12 X(I) = SINF (A(I)) 
Y¥(I) = SQRTF (B(I)) 
Z2(I) = ACI) ** 2 * B(I) ** 3 
GO TO 15 
14% X(I) = SINF (B(I)) 
Y¥(I) = SQRTF (A(I)) 
Z(I) = BCI) ** 2+ A(I) ** 3 
15 CONTINUE 
END 
SUBROUTINE TEST (A, B, J) 
IF (A - B) 5, 5, 6 
J=0 
RETURN 
6 J=i1 
END 
‘END 


If you found significant mistakes, review 80 - 84; otherwise, 
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How do you feel after the last example? If you don’t feel so good, remember----- 


Rome was not built in a day, 
Rome was not torn down in a day, 
But Rome was conquered} 


Once more we emphasize the relationship between arguments in a CALL statement and 
corresponding arguments in a SUBROUTINE name statement. For example! 


CALL JOSEY (L(I), X + Y, 2.6) and SUBROUTINE JOSEY (M, A, X) 
(1) There are 3 arguments in each--the NUMBER of arguments must be equivalent. 


(2) L(I) and M are FIXED POINT arguments, X + Y and A are FLOATING POINT, and 
2.6 and X are FLOATING-~CORRESPONDING arguments must be of the same mode. 


(3) L(I) is a FIXED POINT variable, X + Y is a FLOATING EXPRESSION, and 2.6 
is a FLOATING CONSTANT -- the CALL statement arguments can be FIXED or 
FLOATING VARIABLES, CONSTANTS, or EXPRESSIONS, 


(4) THE SUBROUTINE arguments are DUMMY VARIABLES, As such they can only be FIXED 
or FLOATING pop-subscripted VARIABLES, These DUMMY VARIABLES take on the 
values of the CORRESPONDING arguments in the CALL statement. 
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The following three examples summarize the relationship between arguments of 
CALL statements and SUBROUTINE name statements. 


CALL 
(1) 
(2) 


CALL 
(1) 
(2) 


CALL 
(1) 


MATE (X, 5, ZEB) and SUBROUTINE MATE (A,N,C) 
Same number of arguments (3 each). 


Corresponding arguments are of the same mode (X and A, FLOATING; 
5 and N, FIXED; ZEB and C, FLOATING) 


BENTLY (LOT (I), A + B) and SUBROUTINE BENTLY (JOT, VETO) 
Same number of arguments (2 each). 


Corresponding arguments are of the same mode - LOT(I) and JOT 
are FIXED, A + B and VETO are FLOATING. 


FIX and SUBROUTINE FIX 


Number of arguments--NONE, This SUBROUTINE needs no arguments since 
it is constructed to be self contained. It may perform an action 
completely unconnected to the main program, If it must communicate 
in some way with the main program, this can be done through the use 
of COMMON storage which is explained later. 
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If all the "Smiths" in a city lived in one house, two forms of “sharing” would 
be present. First all Smiths share the same house or location; second, all share the 
same name. 

Efficient programming is possible if the same two forms of sharing exist between 
sub-divisions of a program. For example, if programs, JOHN and JANE each use the same 
variables--why not share these variables? Or if JOHN and JANE each use temporary 
locations in arriving at final results--why not share these temporary locations? 


The Fortran COMMON statement provides for sharing variables and locations, as; 


Example 1 


Main program, WALK, and subroutine, 
RUN use the same variables A and B. 


Example 2 


Main program, AWAKE, and subroutine, SLEEP 
use the same temporary array, DATA. 


PROGRAM WALK 
COMMON A, B 


SUBROUTINE RUN (X, Y) 
COMMON A, B 
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PROGRAM WALK 
DIMENSION DATA (1000) 
COMMON DATA 


SUBROUTINE SLEEP (X, Y, D) 
DIMENSION DATA(1000) 
COMMON DATA 
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When COMMON is used in subedivisions of the same program, assignments in COMMON 
storage are made from the same reference point. Thus in the example below, COMMON B, 
in subroutine GENER, assigns B to the same locations as were assigned Y by COMMON Y 
of main program. One can see that relative order of variables and of array names 
is important. Also note that array names are listed as non-subscripted variables. 


PROGRAM PIECE) 
DIMENSION X(200), Y(200) 
COMMON Y 

CALL GENER (X) 


SUBROUTINE GENER (A) 
DIMENSION A(200), B(200) 
COMMON B 

DO 15 I =1, 200 

A(I) = B(I) + ACL) 
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In this example, both the main progran, 
PIECE, and the subroutine, GENER use the 
same 200 Y values. Thus 200 storage 
locations are shared through COMMON rather 
than the main program and subroutine each 
requiring 200 locations. (In the main 
program these variables are called Y and 
in the subroutine they are called B) 


92 
“ELEGANT” use of COMMON can lead to interesting results. For example. 


PROGRAM HUMP PROGRAM MIX 

COMMON X, Y COMMON X, Y 

CALL NOARG CALL MIXER 

END END | 

SUBROUTINE NOARG SUBROUTINE MIXER 

COMMON A, B COMMON Y, X 

END END 

END END 

Subroutine variables A and B, are Subroutine variables Y and X, of the 


assigned same locations held by variables subroutine are assigned to respective 

X and Y of the main program. Thus one is locations previously held by X and Y of the 
able to transfer results from a subroutine main program. The result is to interchange 
to the main program without having to list names’ of variables from main program to 
arguments of the subroutine. subroutine. 
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In summary, the COMMON statement is given by simply writing "COMMON" followed 
by a list of FIXED or FLOATING VARIABLES. If the variable name represents an 
array, the array name must appear without a subscript. A DIMENSION statement 
given previously provides for the array range. Thus if the array A(200,200) is 
to be assigned to COMMON, one would write: 


DIMENSION A(200,200) 
COMMON A 


Keep in mind, each COMMON statement operates within its own individual program-- 
that is, a sub-program, such as a SUBROUTINE, will require its own COMMON statement. 
Also, the order and number of assignments is crucial. Since programs and 
sub-programs are compiled independently, the use of more than one COMMON statement 
within the same program or same sub-program will make additional assignments in 
succession without overlapping previous assignments. However, COMMON statements 
appearing in each sub-program will make assignments which start from the same COMMON 
reference point. The degree of overlapping will depend upon the list order of the 
COMMON statements. 
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Try the following exercise! 


There are 3000 values stored at array, TANG. Program, TANGENT, computes the 
tangent of each value and stores the function results, Y, at the same storage locations 
reserved for TANG, The following program was written using the tangent function 
subroutine, TANF, from the Reference Library. 


PROGRAM TANGENT 
DIMENSION TANG(3000), ¥(3000) 
COMMON TANG, Y 
DO 20 I = 1, 3000 
20 Y(I) = TANF(TANG(I)) 
ED 
END 


If you think the above program is correct, Go To 96. 
If you think it is incorrect, Go To 95. 
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You deserve a Gold Medal for excellance! The solution is incorrect, and even 
the irrevelent Y (for results), given in the problem statement, did not cause you 
trouble. The correct solution would not require the use of Y or COMMON, and could 
be written as; 


PROGRAM TANGENT 
DIMENSION TANG(3000) 
DO 20 I = 1, 3000 

20 TANG(I) = TANF(TANG(I)) 
END 
END 


Try writing program, MEASURE, described below. 


There are 1000 positive two digit numbers at array, DEBIT. A subroutine, MEAN, 
computes the average, AVER, of all the numbers in the array that are less than 50. 


The main program perserves DEBIT, and forms a new array, BALANCE, by adding the 
average, AVER, (found in the subroutine) to the first 500 numbers of DEBIT and sub- 
tracting AVER from the last 500 numbers of DEBIT, 
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The stated problem wanted each initial value at TANG to be replaced by a new 
value. In Fortran, a procedure such as this is easily accomplished by writing the 
variable to be replaced on the left side of the equal (=) sign and the modified form 
of the variable on the right side. For example, to replace each Z by Z + A, one 
simply writes, 

Z=Z+A 

Therefore, in the previous example, you should have recognized the irrevelency 

of Yeeknowing that to replace TANG, all that was needed was a statement similar to: 


TANG(I) = TANF(TANG(T)) 


The second misconception may have come from a misunderstanding of COMMON. The 
statement "COMMON TANG, Y* does not assign arrays TANG and Y¥ to the same location. | 
Rather, it assigns EACH to locations--in the sequential order: array TANG, array Y. 
These arrays may then be referenced by other COMMON statements that appear in 
sub-programs (subroutines) of the total program, In the given example, there were 
no sub-programs,and COMMON had no useful purpose for being there. 

Now, try re-writing the previous example, correctly and then, 
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PROGRAM MEASURE 
DIMENSION DEBIT(1000), BALANCE(1000) 
COMMON DEBIT 
CALL MEAN(AVER) 
DO 5 I = 1,500 
5 BALANCE(I) = DEBIT(I) + AVER 
DO 10 I = 501, 1000 
10 BALANCE(I) = DEBIT(I) - AVER 
END 


SUBROUTINE MEAN (AVER) 
DIMENSION DEBIT(1000) 


COMMON DEBIT 
CTR = 0. 
AVER = 0. 


DO 15 I = 1, 1000 

IF (DEBIT(I) - 50.) 6, 15, 15 
6 CTR = CTR + 1, 

AVER = AVER + DEBIT (I) 


15 CONTINUE 
AVER = AVER/CTR 
END 
END 


Note: By assigning the array, DEBIT, to COMMON, 1000 locations are saved since both 
the main program and the subroutine are able to use the same locations. 
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How about another example? 
Assume there are 1000 positive two digit numbers (all integers) at array, LIST. 
A subroutine DIGIT, consecutively adds 10 to each LIST term until the accumlated sum 


LSUM, for that term is greater than 100. The number of additions of 10 for that term 
are also kept in counter, KTR. 


The main program uses LSUM and KTR of each term (found in the subroutine) to 
compute the following for each term: 
(1) The right most digit (units’ digit) of the LIST tern. 


(2) The quotient, J, found by dividing the right most digit(found in step 1) 
by KTR. 


(3) The value, Y, formed by taking the sine of the quotient found in step 2. 


Try writing this program in two ways. first, by making use of COMMON; second, 
by not using COMMON, 
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(using COMMON) 


PROGRAM FACTOR 
DIMENSION LIST (1000), Y¥(1000) 


COMMON LIST, LSUM, KTR Two points are worthy of note! 

DO 50 I =i, 1000 
CALL DIGIT(TS Variable, J, is made, FJ, since 
FJ = (LSUM - 100)/KTR the NEXT STEP involves the SINE 

50 Y(I) = SINF(FJ) of FJ and the sine routine requires a 
END FLOATING POINT argument. 
SUBRUUTINE DIGIT(J) 
DIMENSION LIST (1000) Transfer of data between main 
COMMON LIST, LSUM, KTR program and subroutine is made by 
LSUM = LIST(J) putting LIST, LSUM, and KTR in 
KTR = 0 COMMON and by setting up a needed 

2 LSUM = LSUM + 10 link between I and J in the CALL 
KTR = KTR + 1 and SUBROUTINE statements. 
IF(LSUM - 100) 2, 2, 4 

& END 
END 
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(without using COMMON) 


PROGRAM FACTOR 
DIMENSION LIST(1000), ¥(1000) 
DO 50 I = 1,1000 


CALL DIGIT (LIST(I), LSUM, KTR) 


FJ = (LSUM - 100)/KTR 

Y¥(I) = SINF(FJ) 

END 

SUBROUTINE DIGIT (J, LSUM, KTR) 
LSUM = J 

LSUM = LSUM + 10 


KTR = KTR + 1 

IF (LSUM - 100) 2, 2, 4 
END 

END 


In this method, the communication links 
between main program and subroutine are 
established by the corresponding variables 
listed in the CALL and SUBROUTINE 
statements. Thus LIST(I) goes to J,LSUM to 
LSUM, and KTR to KTR. 


In this example only 2 more locations 
are required by not using COMMON, 
Recalling a orevious example, where use 
of COMMON saved 1000 locations emphasizes 
that COMMON is more useful in some 
problems than others. 
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Handy Rules of Thumb: 

First draw the curves -- then plot the data. 

In debugging any type of program, no corrections can be made 
correctly after | p.m. Friday. 

The correct ccrrections will be self-evident at 9:01 Monday morning. 
"Bugs" which would normally take ome day to find will take five 
when you are in a hurry. 

Experience is directly proportional to computer tims wasted. 
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What’s in a name? 
A rose by any other name 
Would smell as sweet! 


Such said the poet long ago. How about the programmer using Fortran? 


An identification statement is: SUBROUTINE JESSIE (NO, NEVER, MAYBE). 
Can we call “JESSIE” as follows, 


CALL JESSIE (YES, PERHAPS, MAYBE): 


If you think this is 0. K., Otherwise 
Go To 102 Go To 103 
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Men Claim: 


“When a lady says NO, she means yes, When a lady says never, she means perhaps," 
but in this case NO can not mean YES and Never can not mean Perhaps, since FIXED POINT 
and FLOATING POINT variables can not be mixed in the calling sequence of subroutines. 
Arguments of subroutines must agree in mode and number with the dummy variables 
of the subroutine identification statement. 
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Self-evident-appearing facts are often misleading. Thus a lady who puts her 
arm out the car window may not be indicating a left turn--perhaps it only means 
the window is down. The fact that you correctly answered the previous question 
doesn’t mean you know all about SUBROUTINES--so let’s try another! 


Would you say there are one or two rules that appear to be broken in the following? 


Part of Part of 
Main Program | Sub Program 
DIMENSION A(100), B(100) SUBROUTINE EVIDENT (ALPHA, BETA, GAMMA) 


COMMON J, B, K + M 
CALL EVIDENT (X, Y, B, 3) 


If you think there was one rule broken, Go To 105 
If you think there were two rules broken, Go To 104 
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Two rules were broken--the number of arguments must agree, and an expression 
(K + M) can not appear in a COMMON statement. 


"Can I take your order?” asked the counter-man, 

"Yes," I said, "make it a cheese:sandwich and a piece of pie." 

“Pie and cheese sandwich", shouted the man, 

"No," I retorted, "it’s cheese sandwich, then pie--I don’t want pie first". 


Order is important in Fortran too. For example, the list of variables ina 
CALL statement establishes an order of modes (FIXED or FLOATING) which must be 
maintained by the corresponding dummy variables listed in the subroutine identification 
statement. 


The same type of order by modes should be maintained between variables of COMMON - 
statements in a main program and its sub-programs COMMON statements. 


CALL JOE (I, J, A) COMMON M, J, E (main program) 
SUBROUTINE JOE (N, K, D) COMMON N, K, G (subroutine) 


Thus, in each set above, the first two corresponding variables are FIXED POINT, 
the third is FLOATING. 
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Perhaps she was making a right turn! There were two rules broken. 

Expressions are not permitted in COMMON statements--only variable names. 
Thus COMMON J, B, K + Mis not allowed. 

The other error occurred in the subroutine calling sequence. CALL supplies 
four arguments (X, Y, B, 3) while the subroutine EVIDENT is prepared to accept 
only three (ALPHA, BETA, GAMMA). Also don’t forget that there must be an agreement 
of mode, thus statements similar to the following would correct this error. 


CALL EVIDENT (XZ, Y, B, 3) 
SUBROUTINE EVIDENT (ALPHA, BETA, GAMMA, N) 
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A main program uses the following statement: 
COMMON A, B, I, D 
A subprogram uses three variables, BETA, GAMMA, and X, BETA uses the data 
assigned to "A" in the main program. GAMMA uses the data assigned as "D® in 
the main program, In addition, the variable "X* in the subprogram is referenced 
through COMMON by the main program. 


Assume the COMMON statement of the subprogram is written as either 1,2, or 3; 


1. COMMON BETA, X, I, GAMMA 
2. COMMON BETA, GAMMA, X 
3. COMMON BETA, B, X, GAMMA 


Choose the correct one above, and depending on your choice (J=1, 2, or 3); 


Go To (107, 109, 108), J 


THE LAW OF THE TOO SOLID GOOF 
In any program, the part that is most obviously correct beyond 
all need of checking is the part that is wrong. 
Corollary I. No one whom you ask will see it either 
Corollary II Everyone who stops by with unsought advice will 
see it inmediately. 
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Shed a tear for Geary, 

Hang the head for Ned, 

COMMON goofs brought reproofs, 
With censure gone, let’s go on! 


A SUBROUTINE, "TASK" is the last subprogram of a larger program. It uses 
three arrays, AID, APT, ard ALERT, each of 400 terms kept in COMMON storage. 
Statement number 62 in the SUBROUTINE, examines the value of the expression 
formed by the sum of the last terms of arrays AID and APT. If zero, control 
goes to statement 72, if not zero, an exit is made from the SUBROUTINE, 
Statement 72 of the SUBROUTINE examines the value of another expression formed 
by the sum of the first terms of arrays AID and ALERT, If not zero, an exit 
from the SUBROUTINE is made. If zero, control goes to a STOP. 


Write as many Fortran statements as possible from the above description. Then, 
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COMMON BETA, B, X, GAMMA 


The statement above is not correct since the FLOATING POINT variable, X, is 
being assigned to a location defined as FIXED POINT, I, in the main program. 


It is very important that yoy keep the order and mode consistent between 
the common statements of the main program and the common statement of the 
subprogran. 


The correct COMMON statement for the example problem would be: 


COMMON A, B, I, D (Main Program) 
COMMON BETA, X, I, GAMMA (Subprogram) 
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COMMON BETA, GAMMA, X 
The statement above is not correct since the variable GAMMA is assigned the 
same storage location as "BY, instead of “D", in the main program. Also since "X" 


is FLOATING POINT, it can not be assigned to "I" which is FIXED POINT. 


It is very important that you keep the order and mode consistent between 


common statements of the main program and common statements of the subprograms. 
The correct COMMON statement for the example program would be. 


COMMON A, B, I, D (Main program) 
COMMON BETA, X, I, GAMMA (Subprograms) 
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SUBROUTINE TASK | 
DIMENSION AID (400), APT (400), ALERT(400) 
COMMON AID, APT, ALERT 


62 IF (AID(400) + APT(400)) 65, 72, 65 


65 RETURN | 
72 IF (AID(1) + ALERT(1)) 65, 75, 65 
75 STOP 

END 

END 


Some variation of the above solution is possible. For example a second 
RETURN statement could have been used. 

We call your attention again to the fact that a RETURN statement provides 
a subroutine inner exit preceding the normal exit designated by the subroutine END 
statement. Also notice that the second END statement is required to indicate 
to the compiler that this is the last statement of the total program. 
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COMPUTING CENTER HAZARDS 


Necktie caught in magnetic tape handlers during rewind 
Consumptive tendencies induced by air-conditioning 
Sleep-walking maintenance men and programmers 


Being run down by spectators 
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Up to this point, approximately three-fourths of the study course has been 
presented. Only INPUT-OUTPUT remains! The following exercises review information 
given up to this time. 

The cards of a Fortran program (excluding Input-Output statements) fell on the 
floor. After picking them up, the following order was found. Can you unscranble - 
the “pick-up order" to the correct order? 


IF (DATA) 5, 5, 3 
PROGRAM DROPSY 


5 END 
3 DO 10 I = 1, 100 
COMMON A, B, X 
DIMENSION A(100), B(100), X(100, 100) 
END 
DO 10 J = 1, 100 
10 X(I, J) = B(I) + ACJ) 
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PROGRAM DROPSY 
DIMENSION A(100), B(100), X(100, 100) 
COMMON A, B, X 
IF (DATA) 5, 5, 3 

3 DO101=1, 100 
DO 10 J = 1, 100 

10 X(I,J) = B(I) + A(J) 

5 END 
END 


Are you wondering why there is a need for the COMMON statement above? If 
you are, so are we! However, it is left there as an example of superfluous 
statements,‘a "freeloader",--most programs usually have a few! 
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Using the flow chart as a guide, write the required Fortran statements. 
GREATER—> (3) 


s+-+O-f= 3} +0 





No @ 
— he Reduce iavance Ga) A 
GREATER—®(6)—>} 5 fee ? G—>|fevaree L by 1 = ee to00 Tes—» 


The NUMBER GENERATOR subroutine indicated above is3 

















Form 
I=N. 31277 


J = N + 46777 


Form generated 
number 
I-J 





Go To 114 


114 


——— ROTATE 90 DEGREES————_—, 


PROGRAM EXAMPLE 


NN = 0 
CALL GENER (N) 

IF (N) 1, 5, 3 
N= oN 

IF (N- K) 4, 5, 6 
CALL GENER (NN) 

N = N +.NN 

GO TO 5 

N=N.7 
L=L+1 

IF (L - 1000) 2, 7, 2 
STOP 
END 


- SUBROUTINE GENER (N) 


=N+e+i 
N %* 5 
J ** 5 
= N - 31277 
= N + 46777 
=sIeoJ 


za oH Ze =e 
| 
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Matrix “HENRY” is placed in COMMON as follows 


Haya Ha,a Ha,3 COCKER E OO OEE Ha, se 


Hays H2,2 Ha, Cooccsccecocces Ha, 50 
1 4 
Hso,52 Hso,2 Hso,3 ececccconcceoce Hso,s0 


Write a Fortran program which will transpose "HENRY" to "GENTRY" -- also in COMMON, 
"GENTRY" will be? 


Ha oa Haya cdc ccccccccccccccccce Hsooa 


Haya Has Coeeoeeeeesseoeseceoooeos Hsosa 


‘ ‘ ) 
Ha ,s0 Hasso Coeceeseeeoeccesecoseoe Hso»50 


As you can see transposition involves changing rows to columns and columns to rows. 
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PROGRAM TRANS 


DIMENSION HENRY (50, 50), GENTRY (50, 50) 
COMMON HENRY, GENTRY 


D010 I= 1, 50 
DO 10 J = 1, 50 

10 GENTRY (I, J) = HENRY (J, I) 
END 


END 


Statement 10 could have been written as follows: 
10 GENTRY (J, I) = HENRY (I, J) 
Arrays are stored in consecutive locations. In the case of multi-dimension 
arrays, keep in mind, that the first subscript is advanced first while the others 
remain constant, etc. Thus the array A(3, 3) would be stored as follows: 


Aas, Aza, Asa, As2, Aza, As2, Aas, Aza, Ass 
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Write Fortran statements which will provide for the following relationships: 


a) SUM= + (ip + 5)* 
=10 


b) 


c) 


d) 


e) 


f) 


d= 1/2 Jur? .T? , where R = cos (@/2) 


Segment = R?cos ~* R R he. (R-h) J 2Rh — h? 


_ 2 (2d? 2 (2d* . 2 (@) 
Length of are =o [1 + 23) - 5 (24 + 7 . 


Find the angle, ~® between two lines whose direction angles areQdisg 25 
iy 2, and 157” 2, and given bys 
Cos ® = Cos@ 1 Cos Q2z + Cos Bs Cos fp? + Cos 7%, Cos V2; 


2(5a°_x° - 6a* bx* + 8a b? x - 16b*) 
FINT = 35 a Vv ax +b tan x 
a 
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(A) SUM =0. 
DO 10 K = 10, 20 
10 SUM = SUM + (J(K)**2 + 5)**K 


(B) D = SQRT (4.*(COSF(THETA/2,))**2- T ** 2)/2. 
(C) SEGMENT = R**2*ACOSF((R- H)/R) - (R - H)*SQRTF(2.°R*H — H**2) 


(D) T= 2,°D/F 
ARCL = G*(1. + 2./3,*T**2 = 2./5,9T**4 + 2,/7.*T**6) 


(E) THETA = ACOSF(COSF(ALPHA(1))*COSF(ALPHA(2)) + COSF(BETA(1))*COSF(BETA(2)) 
1+. COSF(GAMMA(1))*COSF(GAMMA(2) )) 


(F) FINT = (2,°(5.*(A**3)#(X**3) -6,*(A%*2) *B#(X**2) + 8, °A*(B*92) *X 
1-16. *(B**3) )*SQRTF(A*X+B) *TANF(X)) / (35. *A**4) 


Note continuations (digit 1 in column 6) of (E) and (F). 
If you had several differences from the above -- review cards 2-25 concerning 
formula statements and symbols, otherwise 
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The following questions are related to CARD FORMAT. 


ae 
be 


In what columns must Fortran statements be written? 

Can a Fortran statement start in column 10? 

How is a "Comments" statement indicated? 

How is a "continuation" indicated? 

Can one save cards by putting two or three short statements on a card? 
Can one statement of 670 characters be entered? 

What should Card 1 of the program look like? 

How can columns 73 through 80 be used? 

What should the last card of a program look like? 

How many columns are required for the following: TABLE=(3**2(M+N.6)+K)*4 


Does the order of cards (except for first and last cards) make any difference 
in the program? 


Since a key punch operator will probably punch your cards for you, can 
you suggest a way by which she can distinguish your numeric "0" from the 
letter "O"? (Kathy requested this question) 
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Columns 7 through 72 

Yes-blanks are ignored (but a statement can’t start below column 7) 
#Cc* in column one 

Colunin six is not blank or zero -- it contains a digit 1 to 9 
No-only one statement per card is allowed. 


No-the maximum number of characters per statement is 660.--there is a | 
maximum of 66 statement characters per card and a maximum number of 10 cards 
(Initial statement card plus 9 continuation cards) per statement. 


The first card of a program must be an identifying program name card such as; 
PROGRAM EXAMPLE 


Required Tag—t Program name(Up to 7 characters) 

Columns 73-80 are generally used to serialize the cards, but may be used 
for comments if desired. They are ignored during compilation. 

The last card of a program must be an "END" statement starting in column 7. 
23 columns or columns 7-29 (if no spaces are used between characters) 


Yes-order of the cards is crucial since control is sequential from one 
statement to the next. 


A letter "0" is written ] and a zero-numerié, 0, is written 0. 
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Array, NUMBERS, contains 100 positive FIXED POINT values. Some of these are 
large: values. Some are relatively small, It is desired to find the "factorial" of 
these values. Thus "factorial 5" (often written 5!) is (1)(2)(3)(4)(5) = 120. 

The programmer first decides to compare each term of "NUMBERS® with the constant, 
"KONST", If equal to or less than "KONST", the factorial is found by successive 
miltiplication, as in the above example. If the array term is greater than "KONST", 
the following rule is used. 


(n + 3) log n = n log e + log J2n 
where 7 = 3.14159, e = 2.71828, “logis base 10", and n is the "NUMBERS" tern. 


log n! 


HINT: If logn = J, thenn = 10° 


Write program, TORIAL, which will find and store factorials for the array, 
NUMBERS , . 
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PROGRAM TORIAL 
DIMENSION NUMBERS (100) 
DO 50 I = 1, 100 

IF (NUMBERS(I)) 6, 5, 6 


IF (NUMBERS(I) - KONST) 10, 10, 20 
N = NUMBERS(I) 
N=N.1 

(N = 1) 50, 50, 15 
UMBERS(I) = NUMBERS(I) * N 


FN = 10. ** F 


Some variations are possible. The first IF statement provides for factorial 


Also, formula calculations must be done in floating point since arguments 


for the common log and square root subroutines (from Reference Library), and constants 
tw and e, are all floating point. The statement preceeding 50 then changes the 


floating point back to fixed point. 
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Havana Cuba: (A & P release) 


Fido Castro made public a new computer today 
before a teeming. mass at Della Rosa square. 


Experts say a new saga of computer technology 
began with his announcement. 


Although complete details are hazy, it is 
,reported that Castro has ordered human calculating 
prodigies to be entombed as the inner working 
parts of the machine. 


Castro ended his speech by claiming a real break 


through -- the elimination of programmer (bottle) 
necks! 


123 
The following example terminates this review. We hope you have continually 
reviewed past cards and the regular Fortran programming manual when necessary. 


The 15 Puzzle 


The "15 Puzzle" consists of numbers from 1 to 15 arranged 
in a square array, as FIG. 1. The only rule is: “Any one of the 
numbers to the immediate right, left, top, or bottom of the 
blank space can be moved to the empty space. In FIG. 1, only 
15 or 12 can be moved. The puzzle is: "given some array (for 
example FIG, 2), can the numbers in FIG, 1 be moved (by the rule) 
so as to achieve this new array". 


Mathematically it can be shown that there are 16 factorial 
different arrays possible. Exactly one-half (over 10 trillion) FIG, 2 
arrays are impossible to ever achieve. After this puzzle was 
introduced (about 1870) to the public it became very popular. 
So much time was spent on it, that in France, a law was passed 
prohibiting "puzzle 15" attempts during working hours, Huge 
sums of money were won by clever people who bet others (less 
clever) that certain arrays (which were impossible to achieve) 
could not be achieved. 


Then someone came up with a method which indicated whether or not a specific 
array might be achieved. This method is described on the next card. 
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(1) Starting Array (2) "To Be Achieved Array" (3) Lettered Positions 








A neat trick to determine if a particular array can be achieved is as follows: 


ae 


b. 


Ce 


d. 


Assume individual squares are positioned by letters shown in (3). 


Consider the number (call it "X") in position A, of "the array to be achieved", 
Count how many numbers smaller than "X" are in positions higher lettered than 
A. For (2), this would be zero (count the blank as 16). 


Do this for all positions, and add up the counts thus obtained for positions 
A through P. (Total count for (2) is 1) 


If the blank square is one of B, D, E, G, J, L, M, or 0, add 1 to the sum — 
otherwise leave sum alone, 





If sum is even, the array can be achieved. If sum is odd, it can not. 


Assuming a certain “to be achieved array" "MAYBE" is stored in COMMONS write a 
Fortran program-which will follow steps "a" through "e” above to test if achievement is 


possible. 


If it can be achieved, store a 1 at “ANSWER”, If it can’t, store a 2. 
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PROGRAM PUZZLE pO 30 J = 1, & 
DIMENSION MAYBE(4, 4) DO 20 K =1, & 
COMMON MAYBE DO 20 L=1, & 
KOUNT = O IF (MAYBE(I, J). - MAYBE(K, L)) 20,20,10 
IF (MAYBE(1, 2) - 16) 1, 8, 1 10 KOUNT = KOUNT + 1 
1 IF (MAYBE(1, 4) - 16) 2, 8, 2 20 CONTINUE 
2 IF (MAYBE(2, 1) - 16) 3, 8, 3 30 MAYBE (I, J) = 
3 IF (MAYBE(2, 3) - 16) 4, 8, & N = KOUNT/2 
4 IF (MAYBE(3, 2) - 16) 5, 8, 5 IF (KOUNT - 2°N) 50, 40, 50 
5 IF (MAYBE(3, 4) = 16) 6, 8, 6 4o’ ANSWER = 1. 
6 IF (MAYBE(4, 1) - 16) 7, 8, 7 GO TO 60 
7? (IF (MAYBE(4, 3) - 16) 9, 8 9 50 ANSWER = 
8 KOUNT = KOUNT + 1 60 END 
9 DMOIT=1, 4 END 


Chances are your solution will not be the same as the above, but perhaps it does 
provide an indication of the correctness of your solution. Note how the even or odd 
test was made. The remainder is lost on integer division, and thus, an even multiple 
of 2 test can be used to determine if the number is even, 
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FILTER 


Fortran OUTPU Data 
Paraneters Program Results 


The logical flow of most problems can be shown by the above diagram. Up to 
this point, this course has presented the FILTERING characteristics of Fortran 
Programs. The final section of the course will be devoted to the INPUT - OUTPUT 
characteristics, 


Fortran provides for moving information to or from computer storage by special 
and unique input-output statements. Two types of statements are all that is needed 
for the movement of data! One type indicates the medium or device that is to do the 
moving and at the same time, lists the data to be moved. The other type formats or 
describes (by letter codes) how iow the data is to be moved. These two types may be 
referred to as? 


(1) Input-Output Data List Statements 


(2) Format Statements 
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This Fortran course is oriented, for Input, toward Garde and magnetic > for 
output, toward cards, magnetic tape, and printer. (Other media are not covered 





There are five different “Input-Output Data List Statements", 


"READ" statement --- for INPUT by CARDS 

“PUNCH” statement --- for OUTPUT by CARDS 

“READ INPUT TAPE” statement --- for INPUT by MAGNETIC TAPE 
“WRITE OUTPUT TAPE” statement --- for OUTPUT by MAGNETIC TAPE 
“PRINT” statement --- for OUTPUT by PRINTER 


The same organizational plan is used by all five to indicate what data is to be 
moved, and how it is to be moved. The example below indicates the method used: 


READ 10, A, B, C, D 


"DATA LIST" -- tells what data is 
to be moved--in this case, A, B,~.C, D 


“STATEMENT NUMBER” -- tells where to find a FORMAT 
statement which in turn will describe how to move 
the data. In this example, the FORMAT statement is 
located at statement 10 in the Fortran Program. 


(rReAD” indicates INPUT by CARDS 
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Each “Input-Output Data List" statement (of which there are five in this course) 


consists of three parts: (1) medium name, (2) statement number, (3) data list, 
These are re-emphasized in the following example. 


UT TAPE , ws Cc, D 


medium name = Magnetic 


Tape No. 3 statement number where FORMAT statement is located 


Note that when Magnetic Tape is involved, a tape unit number must be present to 
indicate which tape is used. In writing programs for the Minneapolis 1604 Service 
Bureau, use tape units 2, 3, or 4, since tape unit 1 is reserved for the Reference 


Library. 
Another example follows: 


PRINT 12, XA, XB, XC 
| O's ate list 


Statement number where FORMAT statement is located 


Medium name = printer 
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Of the three parts of the "Input-Output Data List" statement, the third (DATA LIST) 
4nvolves some restrictions. Let’s look at these. 


1. Only variable quantities can appear--constants are not allowed. (See any 
example below. ) 


2. If individual quantities are to be moved, the variable names simply appear 
in the list separated by commas. (See Ex. 1 below) 


3. If an entire array is to be moved, the name of the array can appear in the 
list as if it were a single variable. (See Ex. 2 below 


4, If only special members of an array are to be moved, or if a special order is 
described, a special notation similar to DO loops is used. (See Ex. 3 below) 


Ex. 1 PRINT 16, X, NET, B 
Ex, 2 PUNCH 202, MATRIX 
Ex. 3 WRITE OUTPUT TAPE 4, 75, (DATA(I), I = 2, 60, 2) 


Ex. 3 tells the Compiler to write on Magnetic Tape Unit No. 4, according to 
the FORMAT located at statement number 75. The array "DATA" is to be written 
starting with DATAz, then DATAse, etc. ....2 to DATAso. Also note comma required 
after DATA(I). 
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A subscripted variable in the data list is probably the most difficult part of 
writing "Input-Output Data List Statements". Another example re-emphasizes this aspect. 


+ 


WRITE OUTPUT TAPE 4, 10, ((B(I,J), J = 2, 50, 2), I = 1, 50) 


The above statement provides for data output on Magnetic Tape Unit 4, according 
to FORMAT statement 10. Thinking of the data list as nested DO Loops help. For example. 


((B(I,J), J = 2, 50, 2), IT =1, 50) 
time: Loop-———> 
Outer op —] 


The data would be written as follows 
Bi, 2, Baye ecove Ba 4,80 
B2,2, Ba,» eeecce B2,50 ete. 
Note: An open parenthesis begins the loop and a closed parentheses followed by a comma 


must follow each loop. No comma is needed, however, after the last closed 
parentheses in the list. (The arrows above indicate parentheses and commas. ) 
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We started the Input-Output section by emphasizing that two types of statements 
were involved. The first, which we called "Input-Output Data List Statements" indicate 
the medium to be used and what data is to be moved. The second type statement required 
is the "FORMAT STATEMENT". This statement permits the programmer to tell the Compiler 
how data is to be moved. . Note, the FORMAT statement causes no action in itself. It is 
simply used as a reference for some “Input-Output Data List" statement. It tells how 
the characters (alphabetic or numeric), representing the values of the variables in 


the "Data List", are converted and positioned during Input or Output. 


Every FORMAT statement includes a statement number, followed by the word "FORMAT", 
which in turn its followed by a FORMAT "specification" enclosed in parentheses. For 
example >< 


10 FORMAT (12, £12.4) 
an mae 
tuispeci fication" 


The crux of the problem is to learn how to write the "specification"! 
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One writes format specifications by using codes made up of letters, numerics, and 
symbols. For this course, the following six code forms are recommended. 





nIw = the number, n, of INTEGER FORM FIXED POINT fields of width, w. 
("n® and “w" are decimal .integer's ) 


! 
nEw.d = the number,n, of EXPONENTIAL FORM FLOATING POINT fields of width, w, 
including “d" significant digits to the right of the decimal point. 
(*n", "w", and "d" are decimal integers and "w" includes "d") 


nX = the number, n, of blanks (or spaces) that are to be processed. 
nH = the number, n, of HEADING (or Hollerith) characters to be vrocessed. 
/ = "the slash", This symbol indicates that the next card or line should 


be read or written. 


() = "parentheses", These are used in special ways to indicate continued 
input-output control. 


Study the above forms carefully, then 
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Let’s start with code forms nIw and nEw.d -- which are defined by an example: 


n = number of fields w = width of each field d = number of significant digits 
of width, w, which are to 
the right of decimal point 


READ 10 BETA, 524) “READ” indicates INPUT by CARDS, "10" locates the 
——™ ( FORMAT statement, "BETA" is a 6 member FLOATING POINT 
10 FORMAT (689.3 | 514) array, JET is a 5 member FIXED POINT array. 






5, indicates there are 5 fields or members 

INTEGER, indicates INTEGRAL INPUT - OUTPUT (Input in this case) 
4, indicates each INTEGER is contained on the card in FIELDS of 
4, COLUMNS each (4 columns = width) 


@H3 


n= 6, indicates there are 6 fields or members 
{E = “EXPONENTIAL NOTATION", - data is expressed with a sign (+), the 
decimal point, the significant digits, the letter "E", the sign of 
the exponent (+ or -), and 2 digit exponent. 
w= 9, indicates each number is expressed in a field of 9 columns. 
d = 3, indicates 3 significant digits of the width are to the right of 


the decimal point. 
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The programmer picks the codes to fit his programs--just as a "“saw-bones” picks 


his medicines to fit his patients! For example: 


On line 1, print a sum (less than 1000), now in KTRS 

On line 2, print a product (less than 10,000), now in KTRP 

On lines 3, 4, etc., print 100 terms of array, DATA; four numbers 
per iine, each with six digits of significance 
to the right of the decimal point. 


PRINT 12, KTRS, KTRP, (DATA(J), J = 1, 100) 
12 FORMAT (13,/, I4,/, (4(E12.6, 2X)/)) 
"I* 4s required for KTRS and KTRP since these are FIXED POINT variables. 13 provides 
for an INTEGER field width less than 1000; I4 for a field width less than 10,000. 
Each of the first 2 slashes (/) causes the starting of a new line. 


(4(E12.6, 2X)/) in conjunction with (DATA(J), J = 1, 100) does the following’ 


Prints four(4) FLOATING POINT results per line, in EXPONENTIAL (E) form (for example -- 
+.625732E+O4) each in a field width twelve (12); of which six digits (6) are 
significant and to the right of the decimal point, each followed by two spaces (2X), 
until one hundred (J = 100) terms are printed. 
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We hope the previous example enabled you to deduce the following facts about "I" and "E*, 


1. FIXED POINT variables in the data list require the code form, nIw, in the 
corresponding position of the FORMAT specification. Likewise, FLOATING POINT 
variables in the data list require a code form, nEw.d, in the corresponding position 
of the FORMAT specification. 


2. niw in the FORMAT specification indicates that FIXED POINT INTEGERS (with no decimal 
point expressed) are involved. nEw.d in the FORMAT specification, indicates that 
EXPONENTIAL FORMS of FLOATING POINT numbers are involved. 


- . 1262 Bt t is an example of an nEw.d number, where; 


|: Two digits for exponent of base 10 
@- Symbol for sign (+ or -) of exponent 
Letter "E" meaning “exponent” 
One digit for each significant digit 
The decimal point 
Symbol (+ or -) for sign of number 
Thus - . 1262 E + 03 = 126.2 
Note: To determine the width (w) of an nEw.d form, a total of 6 (see circled) 


positions) mast be added to the number of significant digits desired. Thus if 8 
significant digits (d) are required, the width (w) is 8 + 6 = 14 or E14,8, 
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Data appears on one INPUT card in the following field order. ( A = space) 


Card Card 

Field Real Card Field Real Card 

Number Value Form Number Value Form 
1 5 = A&éad5 5 612.7 = +,6127E + 03 
2 120 = A120 6 =.03172 = «=,3172E - 01 
3 326 = fp26 7 512.1 = +,5121E + 03 
4 64 = AGL 8 70125 = .7125E + 01 


One programmer wrote the following: 
READ 21, IOTA, JET, K, L, (ALPHA(I), I = 1, 4) 
21 FORMAT (I4, 214, I4, 4E10.4) 
Another programmer wrote the following: 


READ 312, NNN, MIN, J, NN, (B(I), I= 1, &) 
312 FORMAT (414, 4E10.4) 
If you think both solutions are correct, go to 138, 


If you think only one is correct, go to 137, 
If you think none are correct, go to bed earlier tonight,and start at 133 tomorrow. 
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Actually the two sets of statements given were both correct. The only 
difference being in the naming of variables and the manner in which the part of the 
FORMAT specification corresponding to the FIXED POINT data was written. 


Thus FORMAT (I4, 214, I4) is identical in meaning to FORMAT (414), The 
first, tells the Compiler to read 1 field of 4 columns, then 2 more fields of 
4 columns each, and finally 1 field of 4 columns -- a total of 16 columns. 
The second, tells the Compiler to read 4 fields of 4% columns each -- again 
a total of 16 columns. 


Either set of variable names used was correct since it makes no difference how 
a variable is named as long as the mode of the name fits the mode of the data 
(FIXED or FLOATING), and the name is not longer than seven characters, 
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The next example uses the specification characters: "X", "H" and "/" (slash). 
Mean, M, and sigma, IGMA, (both integers), and 100 random numbers, are computed. Each 
random number of array, RANDOM, contains FLOATING POINT numbers of eight significant 
digits. Write the OUTPUT statements which will list the output on Magnetic Tape Unit 4 
in a form similar to the following ( A= space) (~= an output digit) 


MEANAANDASIGMA = -.-A4A._. (Line 1) 


4. random numbers per line separated 
by 5 spaces each, (Lines 2,3,4, 2. to 26) 


WRITE OUTPUT TAPE 4, 17, M, IGMA, (RANDOM(J), J = 1, 100) 
17 FORMAT (15HMEANSANDASIGMA= 3, x.0Otié«éd2zy /; (4(E14.8, 5X)) 


t 


Provides for 4 
EXPONENTIAL form numbers 
per line, separated by 

5 spaces each, until 

j = 100 








Write 15 Hollerith characters: «J 
“MEANOANDASIGMA =" 


Produces 3 digits for MEAN 
Provides for 3 spaces 
Produces 2 digits for IGMA 


Starts new line 
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From the previous example, one can see‘that codes: X, H, and / (slash) are 
easily learned. 
X is used to skip or provide spaces within a line or card. 
35X on INPUT means: skip 35 columns. 
35X on OUTPUT means: skip 35 spaces. 


/ (slash) is used to start a new line or card. 
/// on INPUT means: start READING the third card or line. 


//{{ on OUTPUT means: start PRINTING, WRITING, or PUNCHING 
the third line or card. 


Note that a slash (/) performs the same function as a 
carriage return on a typewriter. 


H is used to provide headings (hollerith) characters either on OUTPUT or 
within a FORMAT statement on INPUT. 


7HMINNY*5 on OUTPUT gives: MINNY=5 


7HMINNY=5 on INPUT substitutes the next 7 characters read in during 
INPUT for the 7 characters, MINNY=5, in the FORMAT statement. 
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Joe Skinflint was known as a "tightwad"., One night. when the computer center was 
empty (just before a holiday), he sneaked into the center and ran the program described 
below. | | 


Program, MOOCHER, produces “calling cards” for Joe by punching his name on each 
of 500 cards. Each card will appear as follows, starting at column 34, (4 = space) 


Joe AA Skinflint 





Write the program, MOOCHER, and don’t worry about the punched cards being 
interpreted---that’s Joe’s problem. 
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PROGRAM MOOCHER 
pO 5 I =1, 500 
5 PUNCH 10 
10 FORMAT (33X, 14HJOE 44 SKINFLINT) 


Provides for punching 
14 Hollerith characters: 
JOESA SKINFLINT 


Provides for 33 spaces 
before punching starts in column 34 


Notice that the data list statement, PUNCH 10, must be placed in a DO LOOP 
since the data list has no variables. In other words, one can not write, 
"PUNCH 10, I = 1, 500° since no subscripted variable is possible for "I = 1,500" 
to operate on. 
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The most difficult aspect of writing FORMAT specifications is the correct use of 
parentheses for subscripted variables. We have been using some of these concepts in 
past examples, but more details are now given. 


Three basic "parenthesis structures" are important. 


B : 
Exafiple: (E14.8). Process left to right, to 
end of list. Each time closed parentheses is 


A BS repeat to end of list reached, a new line of output is started. 
eo, 


Examples (13,(E6.2)). Process left to right 
C to inner most parentheses. Remain within the 

inner most parentheses to end of list. Each 

time inner most closed parentheses is reached, 


AO then Bn epee eee a eeeesen ane ttin at te 
to end of list t parentneses,. 


Example: (215,6(E14.8),/). Process from left 
to right to first inner parentheses. Remain 
within the inner parentheses for 6 lines; then 
continue processing to the next closed paren. 

n theses. If data list is not exhausted, start 
A—eB, B—»C "n"™ times, again from first open parentheses on a new line 
C—eD, repeat from A if list of output.. 

is not exhausted. Go To 143 


Wanted: 


(from the Los Angeles Post) 


programmer for new IRC computer 
system. must be experienced with 
other systems. will give the right 
man a free draining course which 
will enable the IRC to be programmed. 
for right man willing to dry up 
salary $25000 after we drain you. 


- i 
Mrs. J. C. Prim supervised the programming section of the Outer Banks project. 
Known as a tireless and skilled FORTRAN "pusher", she insisted upon the utmost 
efficiency from her section programmers (slaves). 


Her annoyance with Herb Waster (who had often previously wasted her time) reached 
a breaking point when she examined the following portion of a program he had written. 


WRITE OUTPUT TAPE 3, 8 (BETA(I), I = 1,5) 
8 FORMAT (E14,8) 
PRINT 22, (ALPHA(I), I= 1,5) 
22 FORMAT (5(E14.8)) 


"Can’t you see", she snapped at Herb, “that both format statements do exactly the 
same thing -- and consequently one can be eliminated entirely from the program", 
“And don’t argue with me", she continued, as Herb started to protest, "or I°l1l have 
you terminated immediately”. 

What do you think? Was Mrs. Prim right? 


If you think she was, Go to 144, 
If you think she was not, Go to 145. 


1h 


Here lies the soul of Janella Prim 

Alas, for her, life became too grin, 

She withstood slander, rumors, and blame 
But Herb’s "double format" ended the flame. 


There was a difference! 


(E14,.8) provides for one number. per line. 


(5(E14.8)) provides for five numbers on one line. 


Go To 145 


Now that Mrs. Prim met her “Wasterloo", let’s try the following: 


An engineer wants to PRINT the OUTPUT of one of his problems in a 
form similar to the following: ( A = space,——= output digits) 


(20 spaces) 
DISTANCEA ...cccYeccccccccee ee & STRESS 
+2 line spaces 


<«— 200 ————— 
ome Results 


TTT 


DISTANCE and STRESS are headings (Hollerith) for computed results. 
The 200 results are each computed to 3 significant figures. 2 line spaces 
exist between headings and first lines of computed results. 


Write OUTPUT statements which will provide for PRINTED pages in a 
form similar to the above, then 
Go To 146 
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The two OUTPUT statements, referenced to the OUTPUT is show below: 


PRINT 63, (DIST(I), STRESS(I), I = 1, 200) 
63 FORMAT (43K, SHDISTANCE, 20X, 6HSTRESS, ///, (42K, 2(89.3, 20X))) 


200 Results 


ves 2 line spaces 
“en ———3 Spaces ——————=g» DISTANCE we—20 Spaces——geST OQ 13 Spaces 


<a————- 42 Spaces ——_——_—» 





200 200 
Results Results 


Gor = 447 


| 


"“Systy" and “Cysty" were identical twins. In physical appearance both were alike 
but in work habits they differed markedly. For example, "Systy" once wrote a program 
using the following statements, 


READ 10, JET, (ALPHA(I), I = 1,100) 


"I’m having trouble", said Systy, “I try to read in the number 123 and I get 
1230, and I can’t compile the FORMAT statement. That computer isn’t working right”. 
"It’s not the computer", said Cysty. “First, 1230 came in because you wrote 
123 as 123A 0n the card. You did not RIGHT ADJUST your input. As for the FORMAT 
error, you forgot a very important rule. The number of open parentheses mst 
equal the number of closed parentheses--you have two opens and only one closed". 
“All I hear are rules", moaned Systy, "don’t you know rules are made to be 
broken”, 





"If that is true", replied Cysty, "then compiler programs are written not 
to run, because the "break a rule Compiler “has yet to be born”, 
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Programs which do not provide for changing INPUT or OUTPUT are restrictive. 
For example, a program written to READ in exactly 100 data cards is not applicable 
for 101 cards, etc. To attain flexibility, one often provides for the first card 
of a "data batch” to be a contro] or header card. This card then regulates control 
of the data that follow. The following example emphasizes this technique. 


An array, FILE, is READ in from cards. Card format is: 


Card 1: Number of cards, N, in the FILE, Assume N can range 
from 1 up to 10,000. 


Cards 2,3,4 --- to N: One number per card ranging from -.999 to +.999 


Input statements which will provide for READING one FILE at a time are; 







READ 20, N, (FILE(I), I = 1, N) 
20 FORMAT (I4,/, (E9.3)) 
eet Me” 






ovides for READING 1 
number per card until N 
cards are READ 


Provides for READING a 
value (up to 4 digits) 
for N, and then startin 
to READ the next card 


Go To 149 
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Format code form, nIlw, means INTEGER FIXED POINT. Therefore, the data list variable 
must be a FIXED POINT variable, it must begin with I, J, K, L, M, or N. 


Format code form, nEw.d would probably be used since most calculations involving pi 
would want to carry some digits to the right of the decimal point in pi. 


The array name would have to be changed to FIXED POINT - for example, JBETA, 
One line of heading, TABLES, will be printed starting at the left margin. 


No. INPUT data cards can use the whole card width of 80 columns. All other INPUT - 
OUTPUT media use 120 characters. 


False. The pattern inferred is 3 fixed point numbers followed by repeated floating 
point numbers until the list is exhausted. The 3 fixed numbers are not repeated. 


-.00612E + 05 = -612. For 6 significant digits, the form would be: §E12.6 
No. All numbers are written in decimal-the compiler does the conversion for you. 
False. 8X provides 8 spaces, spacing to the fifth line is://///. (5 slasnes) 


Yes. ™50X" requires 50 columns, "1266.2" requires 72 columns. The maximum number 
of columns allowed is 120 for magnetic tapes and line printer, and 80 for cards. 


No. . Input data must be RIGHT ADJUSTED within the field. Thus, A234 will equal 
230, not 23. 


SaP=2 (first 5 characters) will replace JONES in the format statement. The 
next 3 characters (177) will be assigned to the variable “JEITY™, 


Go To 151 


151 


Program, MISFIT, was having trouble! Date coming in was incorrect. In checking 
the card format against the Input statements, the following correspondence was found: 


CARD FORMAT 
Column 1 
. . 
ms AAAte — 2 we ee ee ee E 







Exponential Number 


"J" Value Spaces 


CORRESPONDING INPUT STATEMENTS 


35 FORMAT (13, 5X, E16.10) 
Can you see any trouble spots in the above. If not, go to card 153. If you do see 
"trouble spots”, 
Go To 152 
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You were right, there were two trouble spots. First, the FORMAT statement called 
for "I3" (which infers 3 columns for width) and only 2 columns were used--the number 
of the subscript was not RIGHT ADJUSTED. Second, there were 16 columns provided for 
an E field of 20 columns, This means the exponent value would be lost in reading in. 


A 30 by 30 matrix, called K, is ready for output on the PRINTER. However, only 
the principal diagonal members are to be printed in the following form: 


K(1,1) sot oe 
K(2,2) =—— —— etc. 


The output statements follow. DONT'T PEEK--- try writing them first. 
PRINT 27, (N, N, K(N, N), N = 1,30) 
27 FORMAT (2HK(,I2, 1H,, I2, 2H)=, 13) 


Study the above until you understand the solution! Note, the variable in the 
PRINT statement could have been any FIXED POINT variable (for example, "I", "J", etc.) 


Go To 154 
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According to the input statements there were two “trouble spots" on the card. 


(1) The "I" field was not right adjusted--it was not 3 columns wide as 


required in the format statement. Thus the subscript read would be 
10 times larger than intended. 


(2) The "E" field punched on the card was in 20 columns instead of the 16 
columns specified in the format statement and therefore the exponent 
would not have been recognized as part of this number. 


REMEMBER--- 


The field specified in the format must be in exact agreement with the input card 
field, as the example below: 


Format 
Specified Field 


+ Er 


If you now understand the above "trouble spots", you are ready for the 
next problem on 152. 
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At one time or another, we all try to reach “the end of the rainbow". And, 
more often than not, the pot of gold we visualize at its end, seems ever just beyond 
our outstretched fingers. 


Perhaps part of the trouble lies in our inherent belief that all things must 
endeewhen in fact, some may not have an ending. 


Paradoxically, as one appears to approach his various “endings", there is a 
natural desire to look back--to retrospect--and to trace the journey from the 
beginning. Thus the ending flows into the beginning and one finds himself looking 
forward to a new ending. 


In this course of study, we have reached such an ending--perhaps only a 
milestone on the journey toward complete understanding of the concepts already 
presented. Nevertheless, it appears to be a point where one can pause, and looking 
back, ask himself, “have I established a firm foundation of understanding these 
programming concepts?® 


With this question in mind, we hope you will conclude this course by an all 


out attempt to solve the following ten problems. Let them be a barometer by 
which you can judge how close you might be to this course "ending", 
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ONE coming up, NINE to go 
Knock this off and start the show! 
Write a Fortran program for the following hypotheses: 


INPUT -- Cards, three arrays, X, Y, and IOTA, X and Y each contain 1000 items 
ranging from -2.510 through +2.499. IOTA contains 1000 integers ranging 
from 100 to 2000. Each card has the format: 


IOTA(I), 2 spaces, X(I), 2 spaces, Y(I) 
BODY OF PROBLEM --- Form the following AFTER INPUT: 


Cc, =x? + 4] where 
(& 4 4} 4 = 4, 2, 3.022, 1000 


If C, 2.0, Increase C, by IOTA, 


If C, <0, Increase counter, KTR, by 1 
OUTPUT -— By PRINTER. On line 1, (centered), print KTR= and follow this with the 
accumulated count in KTR. Then line space the next 2 lines. Starting 


on the fourth line, print 4 items of C per line, each of 10 significant 
digits, with 5 spaces between items. 


Go To 156 


10 


30 


50 


Notes 


1) The inclusion of array C in the DIMENSION, since 1000 C values were computed. 
2) FIXED POINT IOTA(I) is floated (statement 30) before its use in a FLOATING 
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PROGRAM ONE 
DIMENSION X(1000), ¥(1000), IOTA(1000), C(1000) 
KTR = 0 
READ 10, (IOTA(I), X(I), ¥(I), I’= 1, 1000) 
FORMAT(I4, 2(2X, E10.4)) 
DO 40 I = 1, 1000 
C(I) = X(I) ** 3+ ¥ (I) ** 3 
IF(C(I)) 20, 30,30 
KTR = KTR + 1 
Go To 40 
FIOTA = IOTA(I) 
C(I) = C(I) + FIOTA 
CONTINUE 
PRINT 50, KTR, (C(I), I = 1, 1000) 
FORMAT (56X, 4HKTR=,I4,///,(4(E16.10, 5X))) 
END 
END 


expression. 


3) The statement KTR = 0 could have been eliminated since variables are all set to 


zero initially by the compiler. 
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TWO coming up, EIGHT to go 
Do this too, see confidence grow! 


INPUT -- NONE 


BODY OF PROBLEM -~ Starting with k = 1 and S = 0, and successively increasing k by 1 . 
find successive values of S, (by adding each new S to the preceding 
S value) until S is equivalent to or greater than 0.5772157, where. 


n . 
S = 2 oe log k (natural logarithm) 


OUTPUT -- PRINT out on line 1, starting at the left hand column, the following: 
NUMBER OF TERMS NEEDED = _ 


where ____ contains the last value of k used 


Remember: e* and log x are available on REFERENCE LIBRARY. 
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AB 
re 
z 


0. 


10 + EXPF (.FK) * LOCF( FK) 


-0 5772157) 20, 30,30 
FK + 1, 
1 


ron 
+ 


40 


K 
9 
2S5HNUMBER A OF ATERMSA NEEDED4S =A, 13) 


Notes 


The arguments for Reference Library subroutines EXPF and LOGF must be FLOATING 
POINT; therefore, you may not write EXPF (-K) and LOGF(K). Thus FK is used 
instead of K, 


However, the PRINTED value of K is desired in FIXED POINT; therefore, it is 
necessary to "re-fix" the floated FK so that it returns to FIXED POINT (see 
statement 30). 
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THREE coming up, SEVEN to go 
Don’t give up, know, know, know! 


INPUT -- The corresponding slopes of paired lines are in arrays "S" and "T". The 
number of pairs is 1000, and slopes range from -.99 through | + oI. This 
data is recorded on Tape Unit 3, where each record lists an “s “slope, 

2 spaces, and a “TIT slope. 


BODY OF PROBLEM ..-.- 
_ i = cos O+ 1 + Sj 
X,=2sin 4, ttt 
FIND ¥,=2cos M, where ( & ~ 0,125 4 radians 
25 ~ 218 i Si and Ti are corresponding slopes 


i 4 4, the ANGLE above, is to be written as a SUBROUTINE, using COMMON storage for 
arrays, S and T. 


OUTPUT -- On Tape Unit 4, list on each line, Xi, Yi, and Zi, each to 8 significant 


digits, separated by 5 spaces. Following output, tape 4 is to be rewound 
to its starting point. 
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PROGRAM THREE 

DIMENSION S(1000), T(1000), X(1000), Y¥(1000), 2(1000) 

COMMON S, T 

READ INPUT TAPE 35 10 (S(I), T(I), I = 1, 1000) 
10 FORMAT (E8.2, 2%, 

DO 201 =1, 1000 

CALL ANGLE (I, B) 

X(I) = 2. * SINF(B) 

Y¥(I) = 2, * COSF(B) 
20 2(I) = 2./B 

WRITE OUTPUT TAPE 4, 30, (X(I), ¥(I), 2(I), I = 1, 1000) 
30 FORMAT (3(E14, 8, 5x)) 

REWIND 4 

END 

SUBROUTINE ANGLE (J, BETA) 

DIMENSION S(1000), T(1000) 

COMMON S, T 

FJ = J 

BETA = (COSF(.125 * FJ) +1. + S(J))/(1. + T(J)) 

END 

END 


Note 


(1) The exact correspondence between the modes of the parameters of the CALL statement 
and the dummy variables of the SUBROUTINE statement. 


(2) The argument for the COSF subroutine could have been COSF(THETA) if the statement 
THETA = FJ * 0.125 had been included. ; 
1 
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FOUR coming up, SIX to go 
Hang on now or stub your toe! 


The radius of gyration, R, is given by the rule: 


mass 


2 2 2 
rR? = ™a ra + m2 2 + -<--- By FT, m 


where 





mM, + m2 + ----- Mm, r radius 


Groups of masses and corresponding radii (each of 10 significant digits) are READ 
in as a job group. Preceding each job group is a control card containing the number 
of cards, N, in that group, and the job group number, NJOB. Control card with N= 0, 
indicates that all data has been READ in. Maximum number of cards per job group is 1000 
and NJOB numbers can range from 1 to 9999. 


Write a program which will READ in each job group, calculate R (by above rule), 
and write results on Tape Unit 4, on one line of OUTPUT as follows. (A = space) 


NJOBA =4 ASDAA N4=4 RADIUS A OF 4 GYRATION= 


INPUT cards have the following format: 


Control card: N, 2 spaces, NJOB 
Other cards: mass, 2 spaces, radius 


Go To 162 
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PROGRAM FOUR 

DIMENSION FM(1000), R(1000) 

READ 2, N,NJOB 

FORMAT (I4, 2X, I4) 

IF(N) 5, 200, 5 

READ 7, (FM(I), R(I), I= 1, N) 
FORMAT (£16.10, 2X, E16,10) 

SUMMASS = 0« 

RA = 0. 

dO 10 IT=1, N 

RA = RA + FM(I) * R(I) **2. 

SUMMASS = SUMMASS + FM(I) 

RA = SQRIF (RA/SUMMASS) 

WRITE OUTPUT TAPE 4, 100, NJOB, N, RA 
Go To 1 

FORMAT (7HNJOBA =4,14, 4X, SHNO=4,14, 19HRADIUSA ORAGYRATION=, £16.10) 
END 
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16§ 
FIVE coming up, FIVE to go 
Ask for trouble, you sure get moah! 
Given a 3} x 3 matrix, A,;,. The problem is to develop the derived matrix, Bre, 
using the rules, If row number, r, of Bre is equal to its column number, c, then use: 


k = (r-1) 
Bro = Arc = (Ark) (Ai) ,if k € 1, summation =O 


If row number, r, of Bre is larger than its column number, c, use: 


ss (co1) 
Bre = Are - (Ark) (Ake) ,if k € 1, summation = 0 


If row number, r, of Bre is smaller than its column number, c, use; 


= reo1) 1 if k¢€ i, 
re [or : ~ Or) (he) Gy 7? summation = O 














Arc with nine elements ranging from B(1,14= B(1,2)4= B(1,34= 
00.1 to 99.9, is READ in row by row, where 
one element is on each card. Compute the B(2,14= B(2,2)4= B(2,3)4= 
derived matrix, Bre, by the above rules, 
and PRINT out the derived elements in the B(3,1)A= B(3,2)4= B(3,3.4= 
form: 
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PROGRAM FIVE 
DIMENSION A(3, 3), B(3, 3) 
COMMON A 
READ 10, ((A(I,J), J = 1,3), I = 1,3) 
10 FORMAT (E9.3) 
DO 400 IR = 1, 3 
DO 400 IC = 1, 3 
IF (IR -IC) 100, 200, 300 
CASE 1; ROW IS SMALLER THAN COLUMN 
100 CALL SUMMA (Ik,IC, SUM, IR) 
B(IR, IC) =(A(IR, IC) - ” SUM) /A(IR, IR) 
GO TO 400 
CASE 2: ROW EQUALS COLUMN 
200 CALL SUMMA (IR,IR, SUM, IR) 
GO TO 310 
CASE 32 ROW IS LARGER THAN COLUMN 
300 CALL SUMMA (IC, IC, SUM, IR) 
310 B(IR,IC) = A(IR,IC) - SUM 
400 CONTINUE 
PRINT 500, ((I, Jy B(I,J), J = 1,3), I= 1,3) 
500 hen (3(2HB(,11,1H,11,3H)A =,59.3)) 
E 
SUBROUTINE SUMMA (J,L, SUM, IR) 
DIMENSION A(3, 3) 


COMMON A 
SUM = 0. 
IF (J-1) 10, 10, 20 
10 RETURN 
20 KK = Jel 
DO 30 K = 1, KK 
30 SUM = SUM + A(IR,K) * A(K,L) 
END 
END 
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SIX coming up and FOUR to go 
Last lap boy, row, row, row! 


The formla for linear coefficient of correlation is; 


_ Y (XY) = N(Mx) My) 


J Tex - nao] Ly? - way)? ] 


X:;and Y,are paired members of two groups, X and Y , N is the number of pairs, Mx and My 
are the MEANS of the respective groups. Summations shown are 1 to N. 


Input is from Tape Unit 2, consisting of N paired values of X and Y where each record 
is made up of X, Y; X, Y; X, Y; X, ¥$ X, Y; and X, Y (6 pairs). Each X and Y have 

4 digits of significance. "N" is defined by the initial control card. A maximum 

of 10,000 pairs can be analyzed. 


Output is on Tape Unit 3. Mx, My, and r are listed under the respective headings 
shown below, (Headings shown start at left most column of record) 


CORRELATIONA ANALYSIS 
OAMEAN AX AADAAAAAEANS Y AA AA A ACORRELATION 
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PROGRAM SIX 

DIMENSION X(10000), Y(10000) 

READ INPUT TAPE 2, 10, N, (X(I), ¥(I), I =1, N) 
10 FORMAT (15,/,(6E10.4)) 

SUMK = 0. 

SUMY = 0. 

SUMKY = 0. 

SUMK2 = 0. 


pdoi51=1, N 
SUMK = SUMX + X(I) 
SUMY = SUMY + Y(I) 
SUMXY = SUMKY + X(I) * Y¥(TI) 
SUMX2 = SUMK2 + X(I) ** 2 
15 SUMY2 = SUMY2 + Y(I) ** 2 
FN = N 
FMEANX = SUMX/FN 
FMEANY = SUMY/FN 
= (SUMXY - FN* FMEANX * FMEANY) / (SQRTF((SUMX2 — FN * FMEANX**2)* 
(suMY2 - FN * FMEANY**2))) 
WRITE OUTPUT TAPE 3, 20, FMEANX, FMEANY, R 
20 FORMAT (20HCORRELATIONA ANALYSIS »/,2X, 6HMEANA X, 8X, 
1 GHMEAN.Y, 6X, 11HCORRELATION, /, "3(Bi0. 4, &X)) 


-_ 


ED 
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SEVEN coming up, THREE to go 
Fire up boy, make it glow! 


The “least squares” technique applies to the line, Yy = ax, + b, involves solving 
normal equations. 


Na +bZX, -PY, = 0 N = number of paired points, X; and Yj. 

[axh, + be XxX? -F &Y, = 0 where Summations (> ) shown are from i = 1 
toie=N 

For example; 

GIVEN PAIRED POINTS; One can see that One can solve the NORMAL 


Is = 0.50, Ya = 0.38 N= 3, My = 4, EY; = 3.20 POUNDING using Determinanta. 
Xz = 1.00, Ys = 0.82 =x = 7.50,2%, XY = 6.01 , 
Xy = 2.50, Ys = 2.00 ax + by = ¢ one uses 
dx +tey=f | 
NORMAL EQUATIONS are. \ 3a + 4b - 3.2 = 0 K= -ce + bf 
la + 7.50 - 6.01 = 0 ae ~ bd 
X and Y values are READ in by groups. Each group is preceded by y wma + ec = 


a control card indicating the number, N, of paired points in the group. 
When N = 0, the program is ended. N can vary from 2 to 50, X and Y 
values contain 3 significant digits. Write a program which will PRINT out "a" and “b® 


4n a form similar to; 
| Aw AAAAAB= 
Go To 168 
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PROGRAM SEVEN 
DIMENSION X(50), Y(50) 
N 


IF (N) 2, 6, 2 

READ 3, (X(I), Y(I), I = 1, N) 

FORMAT (259.3) 

SUMX = 0, 

SUMY = 0. 

SUMKY = 0. 

SUMX2 = 0. 

po4iI=1,N 

SUMX = SUMX + X(I) 

SUMY = SUMY + Y(T) 

SUMXY = SUMKY + X(I) * Y(I) 

SUMK2 = SUMX2 + X(I) **2 

FN = N | 

A = (=SUMY * SUMX2 + SUMXY * SUMX)/(FN * SUMX2 — SUMK ** 2) 
B = (.FN * SUMXY + SUMX * SUMY)/(FN * SUMX2 - SUMX ** 2) 
PRINT 5, A, B 

FORMAT (2HA=, 89.3, 5X, 2HB=, 59.3) 

Go To 10 

END 

END 
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EIGHT coming up, TWO to go 
Near the end, nice going Joe! 
Simpson’s rule for approximating a definite integral is; 


Sf £(X)dx = > (Ha) +4 f(a x) + 2f(a + Ax) + 4f(a + Wx)+ ---£(b)) 
For example, using a = 1, b = 2, Ax = 0.25, and the integral Si /1 + 4x dx, gives: 
f(a) = £(1) = fit 411) = ./5 
f(a +Ox) = 4f(1 +.25) = 4£(1.25) = 
2f(a + 24x) = 2f(1+2(.25)) = 2f(1.5) = 
f(a + 3qx) = 4f(1 + 3(.25)) = 4£(1.75) | 
Since a + 4Ax = 1 + 4(.25) = 2 = b} the last term is f(2) = fi + 4(2) =.6 


Note, the last term is reached when(a + nAx) = b and that no number (2 or 4) appears 
in front of the first or last terms. 


2 22 
So, S; Ji + Ux dx a (5 + 4 + or +4 +) = 2.636 approx. 


Card 1 contains the valueAx, ranging from -.25 to +.25. The first two fields 
of Card 2 contain "a" and "b" each of which range from -9.99 to + 9.99. Write a 
program which will use Simpson’s rule to approximate and PRINT the results for the 


integral, b 
cos xX 
f. a 






Go To 170 


10 


12 


13 
16 


14 
15 


20 


170 


PROGRAM EIGHT 
READ 10, DELTAX, A, B 


SuM = COSF (A)/A 

= (B - A)/DELTAX 

c = 4, 
KN=1, J 


Hs 
~8 


"FENS 
— 
On 


- 4.) 13, 12, 13 


BOBO 


= ” SUM + TERM 
IF(TERM-ENDTERM) 14, 15, 14 
CONTINUE 

= (SUM * DELTAX)/3. 
PRINT 20, SUM 
FORMAT (4HSUM=, 59.3) 
END 
END 
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C * COSF(A+FN*DELTAX) /(A+FN*DELTAX) 
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NINE coming up, ONE to go 
Do this and you’re in the dough! 


A monthly pay problem involves the following (1) READ IN two cards for each employee- 
the first is the “employee current pay card", the second is the "employee permanent record" 
card; (2) Make pay computations, up-date and PUNCH a new set of "employee permanent 
record" cards; (3) PRINT a “monthly pay record". Assume the following: 


CURRENT PAY CARDS FORMAT PERMANENT RECORD CARDS FORMAT 
COLUMNS PAY ALLOCATIONS ; 


i 
l 
Name of employee | 


COLUMNS PAY ALLOCATIONS 


Name of employee 

Badge Number 

No. of dependents 

Gross pay to date (XXXXX. XX) 
FICA tax to date (XXX.XX) 
FED, tax to date (XXXX.XxX) 
Net pay to date (XXXXX. XX) 


Badge Number 

Hours worked (XXX) 
Hourly Rate (X.XX) 
Problem ends if col. 38 
is blank (interpreted as 
zero) 





The PRINTED record contains. 


BADGE (line 1) FICA TAX = 39% of first $4800 earned to date 
MONTHLY 4 FICA__ (line 2) FED. TAX-= .18 (ross Pay = $54 (No. of 
MONTHLY AFED.TAX (line 3) dependents. } 


Write the PAY PROBLEM, then 
Go To 172 
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PROGRAM NINE 

READ 10, NBADGE1, NHOURS, RATE, ICODE 

FORMAT (20X, 15, 13, 59. 3, T1) 

IF (ICODE) .20, 50, 20 

READ 25, NBADGE2,NDEPS,GROSS, FICA, FED TAX, PAY 
FORMAT (20H 915, 12, 4E13.7) 
IF (NBADGE1-NBADGE2) 26, 30, 2 26 

STOP 

FNHOURS = NHOURS 

GROSSMP = FNHOURS * RATE 

GROSS = GROSS + GROSSMP 

IF (FICA = 144.) 32, 33, 35 

FICAM = 0.03 * GROSSMP 

IF (FICA + FICAM - 144,) 34, 34, 33 

FICAM = 144, -FICA 

FICA = FICA + FICAM 

REFUND = 0, 

GO TO 36 

REFUND = FICA - 144. 

FICA = 144, 

FICAM = 0. 

FNDEPS = NDEPS 

FEDTAXM = 0.18* (GROSSMP - 54, *FNDEPS) 

FEDTAX = FEDTAX + FEDTAXM 

PAYM = GROSSMP -FICAM - FEDTAXM + REFUND 

PAY = PAY + PAYM 

PUNCH 25, NBADGE2, NDEPS,GROSS,FICA,FEDTAX,PAY 
PRINT 40, NBADGE2, FICAM, FEDTAXM 

FORMAT (SHBADGE, 15, /, 1 2HMONTHLY AFICA, E11.5, 
/,  1S5HMONTHLY 4 FED.TAX,E12.6) 

END 


END 
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From: PLATO, THE MGDERN PROGRAMMER; (GRIN and Company, 1961) 


“If anything can go wrong with a program -- it will." 
"Do not believe in miracles -- rely on them," 


“Smile -- tomorrow will be worse," 


173 
TEN coming up, and then you are done 


ee 
But three cards are needed for this one! 4 
The Gamma function is defined by the improper integral [a =Jo =! 6-* ay 
One technique for evaluating it is the following: (accurate approx. 12 decimal places) 


Y= [.)" = e” where Z is defined by the series: 


Case 1, x>6 z= (& - 2) log x- x +? log 2m + 735 -- aghes + Tpbges - TeBT- 


—_—h___ 1 
* Tié8x* ° ad 


First, form D, where D = xtk, and k is the smallest integer that 
Case 2, $< x6 makes xtk> 6. 
a Second, form E, where E = x(xt1)(xt2) ..0e. (De1) 
Third, find F, eam F = k) (from case 1) 


Fourth, find Xx), using ) =F/E 


Case 3, x¢ = [cx = n | where Mx =, and (1dx) is found by 
, {(a-x) sing Case 1 or Case 2 above. 
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Let us assume that 1000 real arguments (one per card), ranging from -5.1 to 15.9 
are READ. We want to calculate GAMMA for each of these arguments and print results, 
4& per line, separated by 5 spaces each, and accurate to 4 significant figures.Using 
the previous rules, we first flow chart the problem: 









Done for 
1000 _ terms 








Find Gamma, Y, 
using Case 1 







Divide results by 
x(x+1)....(xtkel ) 
(See Case 2) 






Studying the chart one can 
see that finding Gamma (Case 1) 
can be made a SUBROUTINE since 
it is referenced twice above. 
Try writing the program using this 
scheme! 
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PROGRAM TEN 

DIMENSION X(1000), Y¥(1000) 
READ 15, (X(I), I = 1, 1000) 
FORMAT (E9.3) 

DO 30, I = 1, 1000 

M = 0 

FK = 0. 

IF (X(I) - 6.) 7,°5, 6 

IF om -.5) 8, 5, 5 

Mo 

ALPHA = > 14159 * X(I) 

X(I) = - X(I) 

GO TO 2 

FK = FK + 1. 

IF (X(I) + FK - 6.) 5, 

CALL GAMMA (X(I) + FK, 20) 
PRODUCT = 1. 

PRODUCT = PRODUCT *(X(I) + FK - 1.) 
FK = FK - 1. 

IF(FK) 12, 12, 10 

Y¥(I) = Y(I)/PRODUCT 

GO TO 3 

CALL GAMMA (X(T) Y(I)) 

IF (M) 20, 30, 

Y(I) = 3. 343 59/ (XCD) * SINF(ALPHA) ) 
CONTINUE 

PRINT 60, (Y¥(I), I = 1, 1000) 
FORMAT (4.(B10.4, 5X)) 

END 


SUBROUTINE GAMMA (X, Y) 


Z = (X - .5) * LOGF(X) =X + .5 * LOGF(6.28318) 
1+ 1/12") - 1./(360. *X**3) + 1./(1260.* 

2X ** 5) - 1./(1680. *X*e7) + 1./(1188, *K**9 

3) - 691./(360360. *x**11) 

Y = EXPF (2) 
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There is a time for everything! A time to live, a time to die,---a time to start 
and a time to end. You have reached the time to end this course. 


If you have doggedly persevered through this course, we congratulate and offer some 
condolence for having to live with repeated examples, rules, and humor that at times was 
probably rather “corny”. 


We started out with one main objective. to teach a minimum FORTRAN course that 
might enable non-programmers to do most or all of their ow programming when the need 
arises, Again, we reiterate that there has been no attempt to teach every and all 
rules and concepts of Fortran. 


Nevertheless, we are confident that having mastered the "lessons of the cards" you 
can do a creditable job of programming any problem in Fortran. Only.experience through 
practice can prove whether we are right or wrong. 


Why not give it a try? 


Go To Hawaii for a rest! 
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REFERENCE TABLE OF FORTRAN STATEMENTS 


General Form 


A=B 

CALL Name (a1, a2, -- an) 
COMMON A, B, etc. 

CONTINUE 

DIMENSION A(n), B(t), etc. 
DO nI=ma, m2, m3 

END 

FORMAT (specification) 

GO TO n 

GO TO (mi, nz, -- np), 4 
IF (g) ns, na, ns 

PAUSE n 

PRINT n, LIST 

PUNCH n, LIST 

READ n, LIST 

READ INPUT TAPE. i, n, LIST 
RETURN 

REWIND i 

STOP n 


WRITE OUTPUT TAPE 4, n, LIST 


Definition 


Value of B replaces value of A 

Go to subroutine named, args.aa, a2, -- ap, 
Assigns variables listed to COMMON 

Ends DO, if last loop statement is a transfer 
Shows max. range (n, t) of listed variables 
Loop thru n; I = m to ma, ms = increment 
Last statement in subprogram or program 
Specifies how I/O data is moved 

Go to statement number, n 

Go to statement number, n 

Go to nay nay n3 aS: g ¢; =, or >0 

Pause, indicate n 

Printer output using FORMAT (n) and LIST 
Card output using FORMAT (n) and LIST 

Card input using FORMAT (n) and LIST 

Tape i input using FORMAT (n) and LIST 
Subroutine inner exit 

Rewind tape i to start point 

Stop, indicate n (terminates program) 
Identifies subr. and args. a1, a2 <-- 2p 

Tape i output using FORMAT (n} and .LIs? 
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